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Abstract 



If the fundamental scale of gravity is of the order of 1 TeV, black holes might 
be produced at the Large Hadron Collider. This work presents simulations of 
black holes and other exotic models of physics beyond the Standard Model - 
super symmetry, extra dimensional models and string theory. Isolated leptons 
with high transverse momenta can be used to distinguish black holes from 
supersymmetry and models of extra dimensions. Z° bosons and photons with 
high transverse momenta allow the discrimination of black holes and string 
resonances. The analysis of visible and missing energy /moment a and event 
shape variables complement these techniques. 
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Introduction 



In 1964 the theory of quarks was formulated by M. Gell-Mann and G. Zweig 
[20, 21, 22, 23] and in 1979 S. Glashow, A. Salam, and S. Weinberg shared the 
Nobel Prize in Physics for their discovery of the electroweak theory [2 1]. Quan- 
tum electrodynamics (QED), the theory of weak interactions, and Quantum 
chromodynamics (QCD), the theory of the quark interactions form the basis 
of Standard Model (SM) of particle physics. Results from particle collider 
experiments, such as Large Electron Positron LEP (at CERN) and Fermilab 
Tevatron agree with the predictions of the SM to a very high precision. One 
of the major milestones in the experimental verification of the SM was the 
discovery of W ± and Z bosons by the SPS, UAl and UA2 experiments at 
CERN [25, 26, 27, 28, 29]. 

The framework of the SM relies on the existence of an unobserved par- 
ticle known as the Higgs boson, which was postulated in 1964 by R. Brout, 
F. Englert and P. Higgs [30, 31]. In the minimal SM the existence of the 
Higgs particle is required to generate the masses of all other particles. Many 
questions, however, remain unanswered. The SM does not explain why there 
are three quark generations, it does not address charge-parity (CP) violation, 
it does not provide a viable candidate for dark matter, and finally it does 
not incorporate the gravitational force. Various models such as extra dimen- 
sions (EDs) [32] and supersymmetry (SUSY) (see Ref. [33, 34] and references 
therein) have been proposed to overcome the drawbacks of the SM. 

The center-of-mass (CM) energy of LEP was ~ 200 GeV. Tevatron's pp 
beams reach a maximum energy of 1 TeV. The upcoming Large Hadron Col- 
lider (LHC) [18], with 600 million collisions per second and a CM energy of 
14 TeV, is supposed to provide an answer on some, if not all, of the open 
questions of the SM. 

The most viable candidate of new physics beyond the SM is SUSY. SUSY 
provides an explanation for the Higgs mass problem, a candidate for cold 
dark matter, and unification of low energy gauge couplings by introducing 
superpartners to SM fields. The lightest supersymmetric particle (LSP) in 
SUSY is an ideal dark matter candidate: a massive, weakly interacting and 
neutral particle. 

Models with large extra dimensions (LEDs) are an alternative to SUSY. 
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Figure 1: The energy of electron-positron and hadron colliders constructed (filled 
circles and squares) or in the planning stage [12]. 

LEDs provide a solution to the hierarchy problem. If they exist, gravity would 
become strong at the TeV scale. SUSY and LEDs are essential ingredients 
of string theory (ST). If ST happens to be the correct theory of physics at 
high energies, and LEDs do exist, the LHC could even start producing string 
resonances (SRs). 

An interesting consequence of the presence of LEDs would be the produc- 
tion of micro black holes (BHs) in particle colliders [35, 36, 37, 38, 39, 40, 
41, 42, 43, 44, 45, 46, 47, 48, 49, 50, 51, 52, 53, 54, 55, 56, 57, 58, 59, 60, 
61, 62, 63, 64, 65, 66] (for reviews, see Refs. [16, 67, 68, 69, 70]) and Earth's 
atmosphere by cosmic rays [56, 57, 58, 59, 60, 61, 62, 63, 64, 65]. Production 
of Kaluza-Klein (KK) excitations of gravitons and SM [71, 72, 73] particles 
would also be consequences of EDs. 

SUSY, LED models, KK excitations and ST lead to different physical sig- 
natures at the TeV scale. It is thus worthwhile to look into the means of 
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comparing these models [74]. Comparisons of SUSY and universal EDs/little 
Higgs models in colliders have been investigated in Refs. [75, 76, 77, 78, 79, 80]. 
However, the literature lacks comprehensive and quantitative comparisons of 
the various models of physics beyond the TeV scale. 

The aim of this thesis is to show, in a quantitative way, how to distinguish 
BH events at the LHC from SUSY, graviton events and SRs. The original 
work of this thesis is based on Refs. [81, 82, 83]. SUSY, KK excitations and 
SRs can be distinguished from BHs by the use of kinematic and dynamical 
quantities. These methods will be discussed in detail in Chapter 4, where the 
original work of this thesis is included. 

This thesis is organised as follows. The next chapter reviews the most im- 
portant features of the SM and the various models of new physics above the 
TeV scale, known problems with the SM, and how new physics beyond the 
SM is supposed to overcome these problems. The physics of particle collisions 
is discussed in Chapter 2. Simulation techniques are presented in Chapter 3. 
Analysis and results are presented in Chapter 4. Finally, the last chapter pro- 
vides some conclusions and discusses open issues. Notations, units, acronyms 
and symbols are included in the appendices. 
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Chapter 1 



Physics of the Standard Model 

and Beyond 

In this chapter we briefly review the basics of the SM and of models of physics 
beyond the SM. The purpose of this chapter is to introduce concepts and 
notations for our analysis of new physics above the TeV scale of Chapter 4. The 
content of this chapter is by no means exhaustive. For a more comprehensive 
discussion of the SM and physics beyond the SM the reader is referred to 
[84, 85, 86] and references therein. 

1.1 The Standard Model 

The SM of particle physics describes the elementary particles and their inter- 
actions. The known elementary particles can be grouped into two categories: 
matter particles (fermions, half spin) and force carriers (bosons, integer spin). 
The particle content of the SM is shown in the figure below. 

The matter particles can be further grouped into quarks and leptons which 
occur in three generations. The quarks are: up (u) and down (d) in the first 
generation, charm (c) and strange (s) in the second generation and top (t) and 
bottom (b) in the third generation. Quark's charges and masses are listed in 
Table 1.1. The leptons are the electron (e), muon (//), tau (r) and the corre- 
sponding neutrinos (i^), where i denotes e,fi or r. The masses of leptons are 
listed in Table 1.2. The neutrinos are chargeless and have negligible mass. The 
electron and neutrinos are stable. Muons and taus have a lifetime of ~ 10~ 6 
and 10 -15 seconds, respectively. Each matter particle in the SM is also accom- 
panied by its corresponding antiparticle: antiparticles have equal mass and 
opposite charge of the corresponding particles. The bosons include the photon 
(7), the W and Z° weak bosons, and eight massless gluons; these particles 
are the force mediators. The photon is the mediator of the electromagnetic 
force, the W, Z bosons mediate the weak nuclear force and the gluons carry 
the strong nuclear force (See Table 1.3.). 
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Fermions Bosons 




Table 1.1: Charges and masses of the SM quarks [1]. Due to the asymptotic freedom 
structure of QCD, quarks and gluons are not capable of independent existence [2, 
3, 4, 5] (and references therein). Therefore, quark masses are measured indirectly 
through their influence on hadrons. The computational scheme used to calculate the 
quark mass from experimental observations is the source of the uncertainty in their 
mass. 





charge 


mass (GeV) 


u 


2/3 


0.002 to 0.008 


d 


-1/3 


0.005 to 0.0015 


c 


2/3 


3.38 to 3.48 


s 


-1/3 


0.1 to 0.106 


t 


2/3 


170 to 174 


b 


-1/3 


4.1 to 4.5 



Table 1.2: Charges and masses of the SM leptons [lj. 





mass (MeV) 




0.511±0.000000013 




105.6±0.000004 MeV 




1777±0.017 
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Table 1.3: Charges and masses of the SM bosons[l]. 





charge 


mass (GeV) 


force 


7 








electromagnetic 


W ± 


±1 


80.39±0.025 


weak 


Z 





91.18±0.0021 


weak 


9 








strong 


Higgs 





> 114.4 





Gluons bind quarks together and are responsible for the production of 
hadrons, which are divided in baryons (quark triplets) and mesons (quark- 
antiquark pair). Mesons are unstable and quickly decay via the strong force 
or weak force into lighter particles. For example, the ir and the K mesons 
decay through the weak interaction and are relatively long lived with lifetimes 
of the order of 10 -8 seconds. The J/^B meson decays via the strong interaction 
and has a lifetime of ~ 10~ 21 seconds. 

The only massless particles of the SM are the photon and the gluons. 
The masses of other SM particles are generated through the Higgs mechanism 
of electroweak symmetry breaking. This model assumes the existence of a 
scalar field, called the Higgs field. Dedicated searches for the Higgs by CERN 
LEP and Fermilab Tevatron have been conducted in recent years. While the 
Higgs particle remains elusive, these experiments have set limits on its mass. 
According to the most recent data [87] the mass of the Higgs is constrained to 
114.4 < m Higgs < 160 GeV or m Higgs > 170 GeV. 



1.2 Physics beyond the Standard Model 

There are strong reasons to believe that the SM is not a complete theory. 
Firstly, it describes the electromagnetic, strong and weak forces but not grav- 
ity. In natural units, the strength of the gravitational force is determined by 
the Planck mass. Max Planck introduced the Planck mass Mpl = ^Jhc/G^ ~ 
10 19 GeV from three fundamental constants: Newton's gravitational constant 
G4, Planck's constant H and the speed of light c. At energy scales of the 
order of the Planck scale, quantum gravitational effects become important. 
This suggests that the SM may be a low energy approximation of some theory 
which is valid at higher energies. The SM is also plagued by the hierarchy 
problem [88, 89]: the fundamental scale of gravity is ~ 10 16 times higher than 
the electroweak scale. Why is gravity so much weaker than the other forces? 

The hierarchy problem can also be understood by considering corrections 
to the Higgs mass. The Higgs mass receives divergent corrections, e.g. from 
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production of a quark- ant iquark pair. The observed Higgs mass vn mggs is 
related to the bare Higgs mass m bare by [90] 

m Higgs ~ m bare ~ ^ i 1 - 1 ) 

where A is the cutoff introduced to regulate the divergence. Assuming that 
the SM is valid up to Mpl, a natural choice for the cutoff is of the order of the 
Planck scale. This would cause the Higgs mass to diverge, unless a fine tuning 
of the order of 1 out of 10 16 is imposed to remove the quantum corrections in 
Eq. (f.f): 

1 m bare _ m Higgs /-, „x 

A 2 A 2 ' 1 ' ; 

The SM also does not address unification of fundamental interactions. 
Gauge coupling constants determine the strength of an interaction, for ex- 
ample the fine structure constant a determines the strength of the electro- 
magnetic force. Coupling constants are not strictly "constants" since they 
vary as the momentum transfer changes in a scattering process. If the SM is 
indeed embedded in some theory valid at higher energies, then at those higher 
energies all the three forces should have the same strength. The SM coupling 
constants, a i} where i denotes electromagnetic (i=l) , weak (i=2) and strong 
(i=3) forces, satisfy [91] 

1 1 h X 2 

where, ai(X 2 ) are the coupling constants at momentum transfer scale X, 
ai(n 2 ) are the coupling constants at energy scale \x (for example the Z° mass 
scale ~ 92 GeV), and hi are numerical factors which depends on the particle 
content of the SM: 6i=7, &2=-3, b 3 =-7 [13]. It is observed that the coupling 
constants do not meet at any scale when extrapolated to higher energies using 
Eq (1.3) (See left panel of Fig.1.1). 



1.3 Super symmetry 

The most studied model of new physics beyond the SM is SUSY. In SUSY, each 
boson (fermion) has a fermionic (bosonic) superpartner. Quantum corrections 
to the Higgs mass due to boson and fermion loops have opposite signs. If the 
SM particles and their supersymmetric partners are degenerate in mass, their 
contributions to the Higgs mass cancel. Thus a Higgs boson with mass ~ 100 
GeV is possible [92] 1 . The presence of superpartners leads to a modification of 
the hi constants of Eq. (1.3). Therefore, in SUSY the running gauge couplings 

1 Note, however, that SUSY is a broken symmetry at low energy and therefore the de- 
generacy in the particle and sparticle masses is not exact. 
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5 10 15 20 5 10 15 20 

log 10 (m/ GeV) log 10 (M/GeV) 

Figure 1.1: Gauge couplings as a function of the energy in SM (left panel) and 
SUSY (right panel). The relative width of each line reflects current experimental 
uncertainties. Reproduced from [13]. 



when extrapolated to high energies meet at the Grand Unified Theory (or 
GUT) scale 10 16 GeV (see right panel of Fig. 1.1.). 

In SUSY, superpartners have identical properties as their SM counterparts 
(masses, charges and quantum numbers), differing only in their spin. The 
super symmetric partners of the gluons (spin 1) are gluinos (spin |), the su- 
perpartners of quarks (leptons) with spin | are squarks (sleptons) with spin 0, 
the partner of the Higgs (spin 0) is the Higgsino (spin |) and the partners of 
the of W,Z and 7 (spin 1) are gauginos with spin |. The gaugino is a common 
name for the supersymmetric partners of the gauge bosons. There are two 
charginos and four neutralinos that occur as a result of mixing between the 
Higgsinos and gauginos. 

As we have mentioned, the Minimal Supersymmetric extension of the Stan- 
dard Model (MSSM) [93, 94] allows for the unification of electromagnetic, weak 
and strong forces at Mqut ~ 10 16 GeV. Since superpartners of SM particles 
are not observed at low energies, SUSY must be a broken symmetry. The 
SUSY breaking scale is generally assumed to be around 1 TeV. A method of 
SUSY breaking which is mediated by gravitational interactions is supergravity 
(SUGRA). The gravitino is the supersymmetric counterpart of the graviton. 
In its minimal version, the physics of mSUGRA is determined by a point in 
the five-dimensional moduli space with parameters: 

• m , the common mass of scalar particles (squarks and sleptons) at Mqut-, 

• To] / 2 , the common gaugino and Higgsino mass at M GUT , 

• Aq the common trilinear coupling at Mgut, 
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• tan P, the ratio of the vacuum expectation values of the two Higgs fields, 

• /i, the sign of the Higgsino mass parameter. 

In the MSSM, the Higgs has five states denoted by H Q , ho, A , H ± . The masses 
of these states are related to the masses of the massive W ± and Z° bosons by 
the relations [95] 



1 

m H ,h - 2 



™a + m l ± V ( m A + m V) 2 ~ Am M m l cos2 2 P\ > (!- 4 ) 
m H ± = m 2 w + m 2 Ao . (1.5) 



Typical mSUGRA parameters which are relevant for LHC processes (LHC 
points) are given in Table 1.4 [96] and Fig. 1.2. 



Table 1.4: MSSM Parameters for typical LHC points. The first five points refer to 
ATLAS and LM1 refers to CMS (see Chapter 2). The scalar mass and the gaugino 
mass are given in Ge V [6j. 



LHC point 


m 


m 1/2 


A 


tan p 




A 


100 


300 


300 


2.1 


+ 


B 


400 


400 





2 


+ 


C 


400 


400 





10 


+ 


D 


200 


100 





2 




E 


800 


200 





10 


+ 


LM1 


60 


250 





10 


+ 



Neutralino gluino g and squark q masses are determined by mo and 
mi/2 through the relations Xi ~ ^1/2/2, X2 ~ Xi ~ m i/2, 9 ~ 3mi/ 2 and 
m(q) ~ (mo 2 + 6m 2 y 2 ) 1 / 2 [97]. The gluino and the first two generations of 
squarks masses do not depend on Aq. The masses of the bosons depend on 
mo and the masses of the fermions depend on mi and Aq [98]. 

A symmetry of the MSSM is -R-parity [93, 94]. The .R-parity of a particle 
is defined as: 

P R =(-lf B + L + 2s , (1.6) 

where L (B) is the lepton (baryon) number and s is the particle spin. Super- 
partners have Pr = — 1 and SM particles have Pr = 1. As a consequence, 
i?-parity implies that SUSY particles are always pair produced in the decay 
of SM particles. -R-parity conservation ensures that a SUSY process at the 
LHC ends in a state with SM particles and two lightest stable SUSY parti- 
cles (LSPs), which are generally neutralinos. Being colorless and chargeless 
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m (GeV) m (GeV) 

Figure 1.2: Five ATLAS points (large dots) shown in the niQ — min plane along 
with lines of constant squark and gluon mass. The shaded regions are excluded 
because of theoretical/ experimental constraints. Reproduced from [14]- 

the LSPs escape the detector and are the source of missing transverse mo- 
mentum If T an important signature of SUSY events 2 . Throughout this thesis 
we will assume that i?-parity is conserved, in agreement with the mSUGRA 
scenario. The mass spectrum for the typical LHC points of Table 1.4 is given 
in Table 1.5. 

Hadronic SUSY processes are dominated by gg, gq and qq channels. Ta- 
ble 1.6 shows the SUSY cross section for the LHC points of Table 1.4. Fig. 1.3 
shows a typical SUSY decay chain. This example illustrates the three impor- 
tant characteristics of SUSY processes; production of jets, leptons and missing 
energy. The jets are produced at the beginning of the decay chain from the 
hadronization of quarks and gluons. The end phases of the decay are dom- 
inated by leptons produced by the decay of charginos and neutralinos. The 
missing energy is due to the LSPs. 



1.4 Extra-dimensional Models 

Alternatives to SUSY are extra-dimensional models such as LEDs [88, 99, 100], 
warped braneworlds [89, 101] and universal extra dimensions [102, 103]. In 
these scenarios, the fundamental scale of gravity M* is lowered to ~ 1 TeV by 

2 If i?-parity is not conserved, the missing transverse energy is reduced by the LSP decay. 
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Table 1.5: SUSY mass spectrum in GeV. The first two generations of squarks and 
sleptons are degenerate. Masses are computed using ISA JET (ver. 7.75) [7]. 

the presence of n extra spatial dimensions or by a warping of the metric. 

The idea of extra spatial dimensions was first proposed by Kaluza and 
Klein in the 1920's in an attempt to unify gravity and electromagnetism. The 
recent interest in extra-dimensional models arises from the possibility that 
these models may provide a solution to the hierarchy problem. 

The ArkaniHamed-Dimopoulos-Dvali (ADD) [88, 99, 100] and the Randall- 
Sundrum (RS) [89, 101] models assume the existence of 3-branes. A 3-brane is 
a surface with three spatial dimensions where matter fields are restricted. In 
the ADD model, the gravitons can propagate in the extra dimensions which 
are compactified at small scales (bulk). 

The idea of EDs also find support from ST [104, 105, 106], the most com- 
plete proposal for a theory of quantum gravity as of now. According to ST all 
particles are different vibrational modes of a closed or open string living in 10 
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Table 1.6: Cross section in fb for production of SUSY particles for the LHC points 
listed in Table 1.4 [8, 9, 10, 11]. 
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1444 


642765 


8259 
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Figure 1.3: Typical SUSY cascade decay in a hadronic process (pp). 
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or 11 dimensions. Open strings have endpoints on the 3-brane, whereas closed 
strings describing force carriers or other multiplets are free to move into the 
bulk. 



1.4.1 Kaluza-Klein Reduction 

The effect of the presence of EDs can be illustrated with a simple example: 
the KK reduction on a circle of radius R [15, 107]. The Klein-Gordon equation 




4 Dimensional Space 



Figure 1.4: Compactified extra dimension of radius R (dotted line). The solid line 
denotes J^D spacetime. Adapted from [15]. 

for a scalar field i/i of mass on a 5-dimensional spacetime with topology 
M 4 x Si, where M 4 is the four dimensional Minkowski space and Si is a circle 
(see Fig. 1.4), is 

dt 2 



-Vl + ml)ij = 0, (1.7) 



where V 2 , is the Laplacian in 5 dimensions: 

v 2 = <?L + f- + + $L as) 

5 dx\ dx\ dx\ dy 2 ' 

Here Xi (i— 1,2,3) are the usual three spatial coordinates and y is the coordinate 
of Si. i/)(xi,y,t) is periodic under the transformation y — > y + 2nR. The 
solution of Eq. (1.7) is of the form 

^ = ZZ- 00 M^tW P5y , (1.9) 

where ^5 = l/R (1=±1, ±2 . . . ) is the momentum along the extra dimension. 
The exponential part of the solution follows from the periodicity condition on 
y. Substituting Eq. (1.9) into Eq. (1.7) it follows 

00 / F) 2 I 2 \ 



Therefore, the scalar field ip in 4 dimensions is equivalent to an infinite number 
of KK states states with masses 

I 2 

m 2 = ml + —. (1.11) 
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The above equation shows that a particle in n dimensions can be described 
by a collection of 4-dimensional modes. The lightest field ip$ has mass 
(this is the KK zero mode or ground state) and the excited modes are doubly 
degenerate m(l) = m(— I). The mass of the KK excitation for R~ l <C m$ is, 
in natural units, 

/ 

m ~ m$ + — . (1.12) 
ri 

Processes with energy ~ Rr 1 should be able to probe excited KK modes. 
Thus, Rr 1 is the threshold energy for detecting EDs. 

Models of EDs can be classified on the basis of KK excitations. In the 
universal extra dimension model, all particles experience KK excitations. In 
the flat compactification model (e.g. ADD), or in the warped compactification 
model (e.g. RS), there is a restriction on which particles possess KK excita- 
tions. The work in this thesis is concerned with the two latter models. These 
will be discussed in some detail in Sects. 1.4.2 and 1.4.3. 



1.4.2 ADD 

In the ADD model, all gauge fields except gravity are confined to a 3-dimensional 
brane. Their quanta do not show KK excitations because they do not feel the 
effect of the EDs; only gravitons are allowed to propagate in the bulk. Ac- 
cording to the ADD model, gravity is a strong force in the higher-dimensional 
spacetime, but appears weak to a four- dimensional observer due to its "leakage" 
in the EDs (See Fig. 1.5). In this model, n EDs are compactified on a circle of 




Figure 1.5: Gravitational lines of force between the two test masses leaks into the 
EDs. This leads to gravity being weaker on the brane. Adapted from [15] . 

radius R. The relation between the gravitational constants in 4 and n + 4 di- 
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mensions, G 4 and G n + 4 , can be obtained by considering the 4 + n dimensional 
Gauss' Law [100] 



(Net gravitational flux over a closed surface C) 
v ; (1.13) 

= S , (3 +n )G 4+n x Mass contained within C, 

where Sd = 2n d / 2 /T(d/2) is the surface of an unit <i-dimensional sphere. The 
force between two test masses m 4 and m^, separated by a distance r is 

F A+n = G A+n ^^ for r«R. (1.14) 

For r ^> R, the test masses do not experience the effect of the EDs: the gravi- 
tational lines of force are constrained within the volume of the EDs. Therefore, 
the familiar 1/r 2 force law is obtained: 

F A = G A -^-. 1.15 

R n r 2 V ) 

The relation between the four-dimensional and the higher dimensional New- 
ton's constants, using Eqs. (1.14) and (1.15), is 

G 4+ „ = R D ~ 4 G A = R D ~ 4 Mpl, (1.16) 

where D = n + 4 is the total number of spacetime dimensions. Since R 
in Eq. (1.16) has the dimensions of inverse mass, G 4 + n has dimensions of 
(mass) 2_D . The D dimensional Planck's constant M* is defined as 

M* = G]^. (1.17) 

The first equality of Eq. (1.16) can also be obtained through dimensional 
analysis [106]. The four-dimensional gravitational potential & g u) satisfies 

V 2 $ 9(4) = 47rG 4 p m , (1.18) 

where p m is the mass density. In D dimensions, Eq. (1.18) can be rewritten as 

V 2 g(n+4 ) = 4vrG n+4 p m . (1.19) 

The left hand sides of Eq. (1.18) and Eq. (1.19) have identical dimensions. 
Therefore, it follows 

Ga+u j^ n+3 = ^ 4 ^3 ^ ^ 4 +« = G/±R D 4 . (1-20) 

Using the above equation and Eq. (1.16), the relation between the fundamental 
scales in 4 and n + 4 dimensions is 

Mp L ~ V n M: +2 . (1.21) 
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If the volume of the EDs is large in 1 units, can be of the order of one 
TeV. 

Setting M* ~ 1 TeV in Eq. (1.21) implies that the radius of the EDs is 
R ~ iQ32/n-i9 me ^ ers n= i i m pii es _R=10 10 km, i.e. a large extra dimension 
with a size of the order of the solar system. The case n = 2 is ruled out from 
astrophysical considerations [108]. Gravity has not been tested to scales of 
~ 10~ 6 mm or ~ 10 11 TeV (for n = 3), and hence one might expect to see 
deviations from Newton's laws at this scale 3 . 

As outlined in Sect. 1.4.1, compactification of EDs causes the appearance 
of towers of KK modes. From Eq. (1.11), compactification of n extra spatial 
dimensions of radius R would result in 

mnQ) = ^y/tl + % + .: + %, (1-22) 

where m n (\) denotes the mass of the n — th KK mode 1. Depending on the 
number of EDs, the graviton modes could be very light. For example, the mass 
of the lowest graviton mode is ~ 25 MeV for 6 EDs . Experimental signatures of 
the ADD model would be production of KK excitations of gravitons and virtual 
gravitons. Since the mass splitting of the KK modes, Am = 1/R ~ l0~ 32 / n 
TeV is extremely small, a large number of KK modes would be produced in 
high energy collisions. 

In the ADD model, gravitons are expected to decay producing jets and 
missing energy (Ifi T ), i.e., a jet or photon recoiling against "nothing". Detec- 
tion of such events would establish the existence of EDs. Searches for graviton 
production have been conducted or are already underway: data from LEP 
and Tevatron have been analyzed to search processes like e + e~ — > r )+$ T and 
e + e~ — > Z+$ T , pp — ► jet+$ T and pp — >• j+$ T [ ]■ The combined LEP 95% 
CL lower bounds on the fundamental Planck mass and the upper bounds on 
the size of EDs are listed in Table 1.7 [110]. 

Table 1.7: LEP lower bounds on M+ and upper bounds on the size of the EDs as 
a function of the number of EDs for the ADD model. 



n 


(TeV) 


R (mm) 




3 


> 1.20 


< 2.6 x 10 


-6 


4 


> 0.94 


; 1.1 x 10- 


-8 


5 


> 0.77 


C 4.1 x 10" 


-10 


6 


> 0.66 


; 4.6 x 10" 
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3 Tests of the Gravitational Inverse-Square Law give bounds on the size of an extra 
dimension: currently it is < 44 /im [109] 
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Planck brane 
y=0 



TeV brane 
y=n 



Figure 1.6: RS model of EDs. The fifth extra dimension is warped leading to the 
weakness of gravity on our brane. SM fields are restricted to the TeV brane. The 
graviton can propagate in the extra dimension. 



The simplest RS model assumes the existence of one extra dimension and two 
3-branes. SM fields are restricted to one of the branes while gravitons are 
allowed to exist also in the bulk. In the RS scenario, the hierarchy problem is 
solved by assuming a single warped extra dimension. The perceived weakness 
of gravity is not due to the presence of a large flat extra-dimensional bulk as 
in the ADD model, but follows from the curvature of the extra dimension. 
RS compactification is on an orbifold which is obtained by applying a circular 
compactification on the extra dimension y, i.e. y — > y + 2irr, where r is the 
compactification radius of the extra dimension y, and identifying y = —y. 
The extra-dimensional coordinate takes values between and nr. Thus, the 
resulting manifold can be seen as a line segment of length nr. 

The branes in the RS model are located at the orbifold fixed points y = 
(Planck brane) and y = nr (TeV brane). SM particles are localized on the 
TeV brane and gravity can propagate in the bulk. The Planck scale is related 
to the fundamental scale of gravity in 5 dimensions by 



where k is of the order of the Planck scale. The hierarchy between the fun- 
damental scales is removed by the presence of the warp factor e~ 2nkr . If kr is 



1.4.3 RS 




(1.23) 
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0.1 0.2 0.3 0.4 0.5 0.6 0.7 0.8 0.9 1 

k/Mp, 



Figure 1.7: Mass of the first graviton excitation mi versus c for kr=10 (solid red 
line), 11 (dashed green line) and 12 (dotted blue line). 



in the range 10-12, TeV scale masses can be generated from the fundamental 
Planck scale ~ 10 16 TeV. A field on the TeV brane with cL IliaSS kx n llcLS cL 
physical mass m n of 

m n = kx n e- kr \ (1.24) 

where x n are the roots of the Bessel function of first order Ji(x n ). The mass 
of the lowest graviton mode in the RS scenario is O(TeV). The strength of the 
graviton-matter coupling is determined by the scale A = (1/ \/8ir)Mpx J e~ k ' Kr . 
In the following sections the ratio c = k/Mpi is chosen to be approximately 
10~ 2 [ : ]. This leads to a graviton mass of the order of ~ 1 TeV, which could 
be probed at the LHC (see Fig 1.7). The mass splitting between the first and 
second excited states of the graviton is Am = rri2 — mi ~ 0(1) TeV. 

Real gravitons in hadronic collisions would be produced via qq — > gG, qg — » 
qG and gg — > gG. The experimental signatures would be multijet or jets+photons 
along with missing Et- Relevant processes for virtual gravitons are // — * 
G* and gg — > G* . Gravitons in the RS scenario behave differently from 
ADD gravitons. Being heavy, they quickly decay into SM particles as in 
G* -> //, gg, 77, Z°Z° and W+W~. 



1.5 Black hole production 

In scenarios with LEDs, pp collisions at the LHC could produce TeV-mass 
BHs 4 . A non-rotating BH is described by the Schwarzschild solution in (n+4) 

4 In this thesis we consider only BHs. For reviews on string balls and branes see [39, 112, 
113, 114, 115] and references therein. 
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dimensions: 
where 



ds' 



-R(r)dt 2 + Rir^dr 2 + r 2 dtt 2 n+2 



R(r) = 1 



R 



BH 



n+l 



The radius of the BH as a function of the BH mass M is 



Rbh(M) 



vr AL 



M 



?r(^) 

. (2 + n) J M 



n + l 



1.25) 
1.26) 

1.27) 



Thus, a fundamental scale-mass BH has a radius of the order of 1 in natural 
units. BH radii for different values of n are shown in Table 1.8 for M=4 TeV 
and M*=l TeV. 

Table 1.8: BH radius as a function of n for a BH mass of M=4 TeV. 



n 


1 


2 


3 


4 


5 


6 


Rbh x 10" 4 fm 


3.64 


2.45 


2.10 


1.98 


1.93 


1.93 



BH production in particle collisions would occur as follows. According to 
Thome's hoop conjecture [116], a BH of mass M is formed when an object is 
compacted in all directions such that 

C < 2irR BH (M), (1.28) 

where C is the circumference of the region where the object is compacted 
into and Rbh is the Schwarzschild radius for a BH of mass M. Assuming no 
gravitational radiation emission the black disk (BD) cross section for head on 
collision is 

&(s;n)=irR% H , (1.29) 

where \/§ is the CM energy of the colliding partons. A more realistic model 
assumes some CM energy being lost as gravitational radiation (see Ref. [117] 
for a more detailed discussion) and non-zero impact parameter. Conservation 
of angular momentum implies that BHs formed with a non-zero impact pa- 
rameter are spinning. They are described by the Kerr solution. A spinning 
BH has a smaller radius than a non-rotating BH of equal mass 

. .„. ^_ _ _ -i 

where J is the angular momentum of the BH. To estimate the gravitational 
energy loss, the colliding particles are described as two Aichelburg-Sexl shock 
waves (see Fig. 1.8) [118]; the overlap of the shock waves (region IV) forms a 



8r(^; 



(2 



n) 



M 
Ml 



n+l 



in 



2fJ 2 



n+l 



4R 2 M 



1.30) 
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trapped-surface (TS) which sets a lower limit to the mass of the BH [119, 120] 5 . 
Depending on the model, BH masses range from 60% (TS model) to 100% (BD 
model) of the total center-of-mass energy. 




Figure 1.8: The u and v lines represent two shock waves moving in opposite di- 
rections, colliding at u=v=0. Region I is before the collision. Regions I, II and III 
are linear. Region IV is highly non-linear and curved, a closed trapped surface is 
formed. The BH is formed in this region. 



The cross section at the parton level for the Schwarzschild BH formation 
process is obtained from Eq. (1.29) as 



a ij ^ BH {s,n) 



1 



8r(^±3)i » 

(2 + ra) 



2_ . .1 



-1 



s 
Ml 



n + l 



(1.31) 



where i,j are the incident partons with CM energy s. The parton cross section 
is independent of the type of incoming partons and depends only on M, M± 
and n. If the energy loss due to graviton emission is neglected, the total cross 
section is obtained by summing the contributions from all possible parton 
pairs: 

o- pp ^ BH (s,n) 

= Y] f dx f fi{x',Q)f j (x/x',Q)a i:j ^ BH (s;n) , ( ^' 32 ' ) 

where s is the squared CM energy of the pp collision, Q is four-momentum 
transfer, fi(x',Q) are the parton distribution functions (PDFs), x' are the 
parton momentum fractions and x=x' i x'j = ~. 

H. Yoshino & Y. Nambu and H. Yoshino & V. S. Rychkov [119, 120] have 
shown that the energy loss due to gravitational radiation is significant for large 



3 For an alternative estimate of the collisional gravitational loss, see Ref. [121]. 



20 



impact parameters and large number of EDs. This reduces the BH mass from 
M = \f§ to M = y(z)V§, where z = b/b max , b max is the maximum impact 
parameter and y(z) ranges between ... 1 [122]. This condition puts a lower 
cutoff on the parton momentum fractions. The lower bound on the parton 
momentum fraction is obtained by the condition x m i n = M^ in /[s y 2 ] where 
M min is the minimum mass for BH formation [122]. The inelastic BH cross 
section is 

o- pp ->BH{s,n) 

= [ 2zdz f dx [ ~T fi( X '^Q)fj( X / X '^Q) Cr ij^BH ■ 

(1.33) 



The total BH cross section (1.33) is shown in Fig. 1.9. The cross section 
increases with energy. For 6 extra spatial dimensions and fundamental scale 
M*= 1 TeV, Upp^BH ~ 6 pb at 14 TeV. This would turn the LHC into a 
BH factory with ~ 2 x 10 5 BHs per year at a luminosity of 10 s3 cm" 2 s~H 6 
The BH cross section decreases with the increase in the number of spacetime 
dimensions and for increasing M*. 




Figure 1.9: Left Panel: BH cross section for different numbers of extra spatial 
dimensions as a function of the CM energy. The fundamental scale M* is assumed 
to be 1 TeV. Right Panel: BH cross section for M+=l TeV (dashed) and M+=2 TeV 
(dotted). 



The distribution of BH masses is shown in Fig. 1.10 for different values of 
M* and n. Higher values of M* lead to more massive BHs. 

TeV BHs decay immediately after formation. The decay starts via loss of 
excess multipole moments (balding phase). In this phase, the BH also loses 
mass by emitting gravitational radiation. In the spindown phase the BH loses 

6 Lepton colliders are expected to produce an even higher number of BHs than hadron 
colliders at equal CM energy because the collision energy is not distributed among the 
hadron constituents. 
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Figure 1.10: Left Panel: BH mass for different values of M+. Right Panel: BH 
mass for different values of extra spatial dimensions, n. Both plots are for 10000 
events. 



angular momentum. This is followed by the evaporation phase where the BH 
loses the bulk of its mass by emitting Hawking radiation [123]. 

The Hawking temperature of a (n + 4) dimensional BH is 

71 + 1 



T, 



4ttR 



BH 



and its lifetime is approximately 



n+3 

1 ( M\ n+1 



T ~ 



(1-34) 



1.35) 



The BH temperature and lifetime for a BH with mass M=4 TeV and M*=l 



Table 1.9: BH temperature and lifetime as a function of n for a BH with mass 
M=4M*=4 TeV. 



n 


1 


2 3 


4 


5 


6 


T H (GeV) 


86 


192 298 


397 


487 


570 


r (xl0~ 26 s) 


1.05 


0.66 0.53 


0.46 


0.42 


0.39 



TeV are shown in Table 1.9. Higher-dimensional BHs are hotter and have 
shorter lifetimes than their four- dimensional counterparts. 

The Hawking process ends when the BH reaches a mass comparable to the 
fundamental scale. At this stage a final n-body decay or remnant production 
(Planck phase) may occur. SM particles produced during the BH decay phase 
are emitted on the brane and can be experimentally detected [124, 125] (see 
Fig. 1.11). As SM particles have overall more degree of freedoms (dofs) than 
bulk particles, brane emission is greater than bulk emission. This result be- 
comes less and less significant in higher dimensions. For example, the power 
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Figure 1.11: Illustration of the Schwarzs child evaporation of a BH, showing brane 
and bulk emission [16]. SM particles are emitted on the brane. The graviton escapes 
both in the brane and bulk. 



loss due to gravitons is negligible in four dimensions whereas it is ~ 25% in 
11 dimensions [126, 127]. 



1.6 String Resonances 



ST is the leading candidate for the unification of the fundamental forces. The 
string scale M s is defined as [106] 

I ^=_L (i.3 6 ) 

where a' is the slope parameter with units of inverse energy squared. The 
slope parameter is related to the four dimensional Newton's constant G 4 [15] 
by 

ot=% (1.37) 
9 s 

The strength of string interactions is controlled by the string coupling g s . 
Combining Eq. (1.36) and Eq. (1.37), one obtains 

M s = g s M*. (1-38) 



Since string effects are expected to appear just before quantum gravity effects 
set in, the string coupling is generally assumed to be of order one. A small 
coupling also justifies the use of perturbative analysis [13]. We have seen that 
in theories with EDs the fundamental quantum gravity scale may be as low as a 
TeV. Thus string effects might appear at the TeV scale. If this is the case, SRs 
would be observed at the LHC before the onset of non-perturbative quantum 
gravity effects such as BH production. String excitation modifications to SM 
amplitudes may even exceed amplitude modifications due to KK gravitons. 



23 



Detection of string events through corrections to SM amplitudes [128, 129, 130] 
at the LHC would be the most direct evidence of ST and EDs. 



In Chapter 4 we will investigate the pp — > string resonance — > 7+ jet process 
[129] 7 . The relevant process for pp —> 7 + jet events is gluon-gluon scattering: 
99 ~^ 91- The string amplitude for this process is [129] 



\M(gg^gj)\ 2 = gtQ z C(N) 



S^SjtjU) Sfi(s,U,t) 



u 



+gtQ 2 C(N){(s^t) + ( 



t 

s «- 



«)}, 



1.39) 



where s, t and u are the Mandelstam variables and 



n(s, t, u) 



u 



v{\-s) r(i-t) 
r(i + t) r(i + s ) 



1.40) 



Here N=3 is the number of D-branes needed to generate the eight gluons of 
the SM, C(N) = jjyj^k is a constant parameter, and Q 2 = |k 2 cos 2 9w ~ 
2.55 x 10" 3 , where k 2 =0.02 and 6 W are the mixing parameter and the Weinberg 
angle, respectively. The values of these parameters are chosen as in Ref. [129]. 



s • 



In the limit s — > nM 2 , the string amplitude possesses poles at n=s/M, 
where n is an integer. The amplitude has the form of Veneziano amplitudes 
[104, 105, 131], a feature common to all string models. The limiting value of 
the amplitude at the n =odd poles is 



\M(gg^ gi )\ 2 = gtQ 2 C(N)- 



t 4 



T(t/M 2 s +n) 
n\) 2 M^[s - nM 2 } \ T(t/M 2 + 1) 



1.41) 



The behavior of the amplitude at the even poles is obtained from Eq. (1.41) 
with the substitutions s — > t and n — + m = t/M 2 in the square bracket term. 
The singularities of the amplitude are smeared with a fixed width T = 0.1 for 
all n > 1 and as 



\M{gg^ 91 )f 



4g*Q 2 C(N) f 

Mf ~ \ (s 
4g 4 Q 2 C{N) 



M 2 ) 2 + (T J=0 M S 
t 4 + u 4 



(s - M 2 ) 2 + (T J=2 M S ) 2 



1.42) 



for n=l [129]. Equation (1.42) includes a correction for spin dependent widths: 
r J=0 = 0.75a s M s and T J=2 = 0.45o; s M s , where a s = g 2 s /AiT is the strong 
coupling constant. The presence of the poles indicates the formation of SR. 
The total cross section for the pp — ► 7 + jet event is obtained by integrating 



7 For a discussion of SR in the dileptonic and diphotonic channels see for e.g. [128, 131, 
132, 133, 134]. 
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the parton cross section over the parton distribution functions of the protons 
[135] 



1 rl 



Tpp^string^+jet = / dx / dx' j dt fx(x', Q)f 2 (x/x\ Q)^, (1-43) 



.r 



where 

da \M(jgg^gj)\' 



;i.44) 



dt 167TS 2 

The limits on t are fixed by the conditions s + u + 1 = and > Pj> min 

where Prmin is the minimum transverse momenta of the two outgoing particles 
of the 2x2 scattering. Figure 1.12 shows the differential cross section of the 




Figure 1.12: Differential cross section of string events for M s = 1 TeV (red dots) 
and M s = 2 TeV (green crosses) with PTmin=50 GeV. SRs are clearly seen when 
s = nMl. 



VV — * 7 + 3 e t process 

^ = y f dx' dt^-h{x/x\Q)f 2 {x\Q)^ (1.45) 

where we have used the CTEQ6D PDFs. The left panel of Fig. 1.13 shows the 
total cross section as a function of Prmin- The string cross section for M s =l 
TeV (solid red line) and the cross section for M s =2 TeV (dashed green line) 
are ~ 5xl0 4 and 10 3 times less than SM, respectively. 

The string cross-sections for F Tmiri = 50 GeV (solid red line) and 300 GeV 
(dashed green line) are shown in the right panel of Fig. 1.13, respectively. For 
Prmin = 300 GeV, the string cross section is ~ 100 times less than Prmin 
= 50 GeV when M+=7 TeV. The choice P<7>mn = 50 GeV leads to a signal- 
to-background ratio of ~ 73. The string cross section is highly suppressed 
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Figure 1.13: Left panel: String cross section for M s = 1 TeV (red dots) and M s = 2 
TeV (green crosses). The cross section for M s = 1 TeV is ~ 44 times larger than the 
cross section for M s = 2 TeV. Right panel: String cross section for the two extreme 
values of PTmin, 50 GeV (solid red line) and 300 GeV (dashed green line) as a 
function of M+ . 



w.r.t. the SM cross section for lower values of 

Prmin'i for example si "" 9 ~ 
10~ 5 for PTmin — 10 GeV. Thus, the discrimination of string events from 
the SM background is more difficult for events with lower PTmin- Both the 
SM background and the signal are substantially reduced for higher values of 
Prmin- For example, at 300 GeV they are reduced by a factor of ~ 98% and 
~ 42% w.v.t. values at PTmin 

= 50 GeV, respectively. The optimal signal-to- 
background ratio is obtained for PTmin ^ 100 GeV. 
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Chapter 2 



Physics of Particle Collisions 



In this chapter we briefly discuss the main features of CERN's LHC and intro- 
duce kinematical and dynamical quantities which will be used in our analysis. 

2.1 The Large Hadron Collider 




Figure 2.1: Illustration of the LHC ring structure showing four of the six exper- 
iments. The energy of the protons is increased in succession by the use of various 
smaller accelerators. Protons, first injected into LIN AC, are emitted with an energy 
of 50 GeV. Next, the PS booster accelerates the protons to 1.4 GeV before it enters 
the Proton Synchrotron (PS). 25 GeV protons from PS is fed to the SPS (Super 
Proton Synchrotron), accelerating them more than 10 times their present energy to 
450 GeV. In the final phase, the protons enter the LHC where they reach their final 
energy of 7 TeV per beam. Adapted from [17] . 

The LHC is a circular proton-proton collider with a circumference of 27 
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km. The LHC is host to various experiments: A Large Ion Collider Experi- 
ment (ALICE), A Toroidal LHC Apparatus (ATLAS) [136], Compact Muon 
Solenoid (CMS) [137], the Large Hadron Collider beauty (LHCb) experiment, 
the Large Hadron Collider forward (LHCf) experiment and the TOTal Elas- 
tic and diffractive cross section Measurement (TOTEM) experiment. The 
ATLAS and CMS (see Fig. 2.2) are general-purpose high-luminosity detec- 
tors designed to cover a wide range of physics from the search for the Higgs 
[8, 9, 10, 30, 138, 139, 140, 141] and SUSY [14, 93, 94, 142, 143] to EDs 
[88, 89, 99, 100, 101, 102, 103]. The main feature of the ATLAS detector is 
its magnet system consisting of superconducting magnet coils. The CMS is 
built with the same goals as ATLAS but it offers an excellent muon tracking 
system. The LHC is based on the barrel plus endcaps design; a cylindrical de- 
tector covers the central region and two flat circular endcaps cover the angles 
close to the beam: ALICE and LHCb detectors are exceptions as they have 
asymmetric shapes. ATLAS and CMS detectors are entrusted with the task 
of studying events with large transverse momentum Pt, a signature common 
to SUSY, EDs and string-mediated interactions. 



The primary purpose of the LHC detectors is to measure particle positions, 
charges, velocities, masses and energies. To achieve this goal, the detectors 
are composed of various sub-detectors. The calorimeter is the most important 
component as it measures the energy of the particles. It encloses the tracking 
system that electronically records the path of charged particles as they pass 
through matter. Calorimeters work on the following principle: they measure 
the ionization energy deposited in the shower when a particle interacts with a 
dense medium (e.g. lead). The electromagnetic calorimeter (ECAL) absorbs 
electrons and photons, the hadronic calorimeter (HCAL) records the energy of 
hadrons (e.g. pions and jets). Photons and neutrons are detected only by the 
calorimeter. Muons are identified by the long tracks in the calorimeter and the 
bulk of their energy is deposited in the muon chamber outside the calorimeter 
(See Fig. 2.3.). At the LHC, the total Pt of the colliding particles is assumed 
to be zero. Neutrinos are weakly interacting particles and cannot be directly 
detected, but their presence can be inferred due to missing energy. 



The LHC follows a right-handed coordinate system where the beam direc- 
tion is along the z- axis and the vertical axis is y. The coordinates can also be 
expressed in terms of (r, 77, 0), where <fi = tan _1 (y/x) and the pseudorapidity 
is r] = — ln[tan(6 l /2)]. The true rapidity of a particle defined as a.ta.iah(p z /E), 
where p z is the momentum along the z direction and E is the energy of the par- 
ticle. The rapidity for massless particles is the same as their pseudorapidity. 
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Figure 2.2: Cross section of the ATLAS and CMS [18] detectors at CERN. 
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Figure 2.3: Passage of charged particles through matter in particle colliders. 
Charged particles leave a track in the tracking chamber, neutral particles shower 
without leaving a track. Neutrinos are undetected. Adapted from [17]. 

2.2 Luminosity 

The performance of the LHC is described by its instantaneous luminosity L, 
which is a function of the number of particles in each bunch(nb), the bunch 
crossing frequency / and the bunch overlap A [1]: 



/ 



A ' 



(2.1) 



The proton beam at the LHC consists of ~ 3000 bunches, each containing 
about 10 11 protons. Particle bunches are squeezed at the interaction points to 
increase the chances of collision to a cross-sectional area of about 10 sq. mm 
[17]. The luminosity can be enhanced by increasing the number of particles 
in each bunch and by increasing the frequency of collision. The latter can be 
achieved by properly focusing the beam at the interaction points. 

The instantaneous luminosity has units of cm~ 2 s~ l . The number of events 
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N for a given process in the time At is 

/■At 

N = a L(z)dz, (2.2) 

where a is the total cross section. The integral in Eq. (2.2) is known as the 
integrated luminosity. The design luminosity of the LHC is 10 34 cm _2 s _1 , to be 
achieved within three years of running. The LHC initial or startup luminosity 
is 10 29 cm" 2 s _1 . 



2.3 Event Shape variables 



Event shape variables are commonly used in particle physics to describe the 
spatial distribution of the outgoing particles. They include sphericity, 2 nd 
Fox- Wolfram moment and thrust [144, 145]. 

• Sphericity: 



S = ^(A 2 + A 3 ), (2.3) 
where A 2 and A 3 are the two eigenvalues of the sphericity tensor 

S «f3 

S<# = ^ P ? ' P \ (2.4) 
s i I Pi I 

Here a, ft = 1, 2, 3 and Pj's are the final state momenta of particles. 
An isotropic event is characterized by values of S ~ 1 whereas S ~ 
for "jetty" events. 

2 nd Fox- Wolfram moment: 



H 2 = Sy^Spi P 2 (cos%), (2.5) 

vis 

where cos#y is the angle between hadrons i and j, E vis is the total 
visible energy of the event and P2{x) is the Legendre polynomial 
of order 2. Values of H 2 closer to zero indicate spherical events 
whereas H 2 ~ 1 indicates that the event is dominated by jets. 

• Thrust: 



T = max| n | = i S ^ n (2.6) 
^i\Pi\ 

where n is the unit vector along the thrust axis. The thrust axis is 
defined as the axis which maximizes projected momenta. Isotropic 
events have T ~ | and dijets have T ~ 1. 
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2.4 Invariant Mass and Jet Masses 



Another two useful quantities in high energy physics are the invariant mass of 
outgoing particles and jet mass. They are defined as 

• Invariant Mass: 

For a system of m particles, the invariant mass is defined as 

M m = (£ E m ) 2 -QTp J 2 . (2.7) 

m m 

where E(p) is the energy (momenta) of each particle. In a 2 x 2 
collision the invariant mass of the outgoing particles is 

M 12 = (E 1 + E 2 f - ( Pl + p 2 ) 2 , (2.8) 
= Ml + Mj + 2{E l E 2 - Pl -p 2 ). (2.9) 

in natural units. Neglecting the masses of the outgoing particles, 
as in often the case in high energy collisions, the invariant mass of 
the system is 

M 12 = y/(E 1 + E 2 f - ( Pl + p 2 ) 2 = yj2 Vm {\ - COS0) , (2.10) 

where 9 is the angle between the two outgoing particles. The invari- 
ant mass is an important quantity in high energy physics because 
it predicts the mass of the particle formed in resonance. 

• Jet Masses: 

The particles of an event are randomly assigned into two groups 
and their squared invariant mass, M\ and M 2 , is calculated. The 
heavy and light jet masses are defined as the larger and smaller of 
Mi and M 2 that minimize M\ + Mf , respectively. 
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Chapter 3 



Simulation Techniques 



This chapter briefly introduces the simulation programs which will be used for 
our analysis: ISAJET (ver. 7.75) [7], PYTHIA (ver. 6.406) [144, 145], CATFISH 
[146, 147] and the string event generator. SUSY simulations will be carried out 
using a combination of ISAJET and PYTHIA. Simulations of graviton processes 
will use PYTHIA. BH simulations will be carried out using the CATFISH Monte 
Carlo (MC) generator. SR and SM simulations will be done with PYTHIA. The 
setup for each simulation is: 

• SUSY: 

— The MSSM mass spectrum is generated with ISAJET; 

— The mass spectrum in SLHA format is fed into PYTHIA; 

— SUSY processes except SM Higgs production are simulated; 

— Unstable SM particles and sparticles are hadronized or decayed with 
PYTHIA. 

• Graviton decay: 

— ff —> G and gg — > G processes are simulated with PYTHIA as 
internal processes. 

• BHs: 

— The cross section of the BH event is calculated in the CM frame; 

— The initial BH mass is sampled from the differential cross section; 

— The BH is decayed through Hawking mechanism and final n-body 
event (or remnant); 

— Unstable quanta are hadronized with PYTHIA. 

• String resonances: 

— The differential cross section for the pp — > 7 + jet process [129] is 
calculated; 
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— The total cross section is obtained from Eq. (1.43); 

— The resonance decay is simulated with PYTHIA. 

• SM background: 

— it, Z°, W, WZ°, WW and Z°Z° SM processes are simulated with 
PYTHIA. 

PYTHIA 

PYTHIA is a FORTRAN 1 MC code used for high energy particle collisions. Typ- 
ically, a high energy collision is a complex process producing a plethora of 
particles. The main purpose of PYTHIA is to generate a detailed event history, 
for example the 2x2 hard-process, initial- and final-state radiation, beam rem- 
nants, fragmentation, decays etc., for a generic detector. Events are generated 
using the following three main steps (see also Fig. 3.1): 

• Hard interaction. The primary hard interaction is generated for process 
under study through the appropriate cross section. The hadrons are de- 
scribed by a set of parton distribution functions f(x, Q 2 ) which give the 
probability of finding a parton with a certain momentum fraction x of 
the total energy of the proton at the energy scale Q 2 . PYTHIA currently 
contains about 300 hard processes which are classified according to the 
number of final particles (2 — > 1 and 2^2 hard interactions) or accord- 
ing to the physical processes. When classifying hard processes based on 
their physical properties, PYTHIA currently supports about 17 different 
scenarios. Some of the processes relevant for this thesis include hard 
QCD processes, W/Z production, new gauge boson production, SUSY 
and extra dimensions. 

• Initial and final state radiation. Incoming- and outgoing- state radiate 
before and after the hard interaction through soft or semi-hard processes. 
For example an initial quark produces a photon as in q — > q^y. 

• Quark- gluon confinement. Quarks and gluons produced in the interac- 
tion are hadronized to produce stable particles. 

Each event in PYTHIA is stored as an event record in a common block called 
PYJETS. Origin, present status, momentum, energy, mass and production ver- 
tices of the particle are recorded. The PYLIST command provides a summary of 
the information stored in PYJETS as shown in Table 3.1. Lines 1. . . 9, show the 
event history, denoted by the exclamation marks around the particle names. 
The first two lines of the output show the initial 7 TeV protons along the z 
direction. Lines 3,4 (5,6) show the initial (final) state radiation particles. Line 

1 Starting with version 8, PYTHIA is written in C++. 
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Figure 3.1: 2x2 interaction showing initial and final state radiation, multiple 
interactions and beam remnants. Adapted from [19]. 

7 contains the graviton produced from gg fusion which decays into a qq pairs 
(lines 8 and 9). The particles produced in the final state are shown from line 
10 onwards. Particle names in brackets signify that these particles further 
decay. The last line of the output shows the total energy and momentum of 
the interaction. 

ISAJET 

ISA JET is another FORTRAN MC for simulating high energy pp collisions. It is 
similar to PYTHIA, but it is primarily used for parton level studies. Events are 
generated as follows in four steps: 

• Computation of the cross section of the hard scattering process and 
simulation of the hard scattering. 

• Simulation of initial state and final state radiation. 

• Hadronization of partons. 

• Jet simulation to describe the "soft" partonic interactions, i.e. multiple 
interactions, beam remnants etc. 
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Table 3.1: Abridged PYTHIA output showing a pp collision to produce a graviton 
(G*). Orig refers to where the current particle originated from. Other symbols have 
the usual meaning. All quantities are in GeV. 



Event listing (summary) 



T 
1 


particle/jet 


orig 


Px 


Py 


Pz 




in 


1 

1 


!p+! 


U 


n nnn 


n nnn 
U.UUU 


Tnnn nnn 
f UUU.UUU 


( uuu.uuu 


n oqo 
U.yoo 


2 


!p+! 


u 


n nnn 

U.UUU 


n nnn 

U.UUU 


^znnn nnn 

-7UUU.UUU 


7UUU.UUU 


n no o 

U.938 


3 


!g! 


1 


1.363 


n Tn o 

U.7U3 


i o/" 1 o r nn 

1868.599 


1 o o o nn 

1868. 6UU 


n nnn 

U.UUU 


4 


!g! 


2 


-1.U25 


U.267 


-27U4.946 


27U4.946 


0.000 


5 


!g! 


3 


2.289 


-1.U8U 


1U97.417 


1U97.419 


0.000 


6 


!g! 


4 


16.514 


13.U34 


-24U.UU5 


240.925 


0.000 


7 


!G*! 





18.8U 


11.95 


857.41 


1338.34 


1027.38 


8 


!u! 


7 


-443.276 


-232.577 


358.721 


615.846 


0.330 


9 


ml 


7 


462. U79 


c\ A A roi 

244.531 


498.691 


h7a c\ A f\f\ 

722.499 


0.330 


10 


(G*) 


7 


18.80 


11.95 


857.41 


1338.34 


1027.38 


11 


(d) 


4 


5.352 


-4.349 


3.821 


7.891 


0.330 


12 


(g) 


4 


1.49U 


-3.33U 


-U.438 


3.675 


0.000 


13 


(g) 


4 


U.253 


U.121 


-U.255 


0.379 


0.000 


14 


(g) 


4 


1.849 


-6.459 


-7.163 


9.820 


0.000 


185 


Pi+ 


137 


-2.183 


U.833 


1.493 


2.776 


0.140 




sum: 




U.UU 


U.UU 


U.UU 


14000 


14000 



ISA JET stores particle information, momentum, origin and decay products, 
in a common block. Information about the production vertex in not stored 
by ISAJET. Contrary to PYTHIA, particles are produced/decayed at the ori- 
gin. Therefore, it lacks detector response simulation. In this thesis ISASUSY, 
included with ISAJET, will be used only to generate SUSY particle masses. 

CATFISH 

CATFISH [146, 147] is FORTRAN MC generator for simulating BH events at the 
LHC. CATFISH interfaces to the PYTHIA MC fragmentation code. Simulations 
using CATFISH are controlled by a set of external parameters and switches: 

• Fundamental Planck scale, 1 < M+ < 14. 

• Number of large EDs, 3 < n < 7. As discussed in Sect. 1.4.2, models with 
one or two large EDs are excluded. BH production in warped scenarios 
such as the RS models (see Sect. 1.4. 3) are not considered in CATFISH 
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since most theoretical studies concerning black holes at colliders have 
been derived for a flat extra-dimensional scenario. 

• Gravitational loss model. Options include the no graviton loss BD model 
and the YN and YR trapped surface models. 

• Minimum BH mass at formation, M m i n in units of M*. 

• Quantum BH mass threshold at evaporation, Q m in ~ 1 TeV. Q m i n is 
defined as the mass of the BH at the end of the evaporation phase. 

• Number of quanta at the end of BH decay, 2 < n p < 18. 

• Momentum transfer model in parton collision. Two choices are allowed: 
M or inverse Schwarzschild radius. 

• Dilepton Invariant Mass Parameters. (See Sect. 4.1.3.) 



CATFISH simulations neglect the energy loss in the balding phase since it 
is poorly understood. Since emissivities of rotating BHs are not known for 
all fields, the description of the evaporation phase is also approximated. For 
simulation purposes, BH generators rely on the Schwarzschild phase. The total 
decay multiplicity in CATFISH is 

_ (n + l)S E» CiKjT^ 

The decay multiplicities per species iVj are 

N l = N-^§^-, (3.2) 

where T-p. and T^. are the relative emissivities of Ref. [126, 127], S is the 
initial entropy of the BH, V s and 1Z S are spin-dependent power and emissivity 



Particle 


Color Flavor Charge Spin 


dof 


Quarks 


3 6 2 2 


72 


Gluons 


8 1 2 


16 


Leptons 


3 2 2 


12 


Neutrinos 


3 2 1 


6 


W,Z 


3 3 


9 


7 


1 2 


2 


Higgs 


1 


1 


Graviton 


1 


1 



Table 3.2: Dof's of the known SM particles, Higgs and the graviton. 
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normalization factors, respectively and q is the number of dof of the particle 
species i. 

Table 3.2 shows the values of c$ for SM particles. The overall dof of a 
species is the product of the particle's color, flavor, charge and spin [1 18]. For 
example, quarks come in 3 color states, 6 flavor states, 2 charge states and two 
spin states. Therefore, the degree of freedom is 3 x 6 x 2 x 2 = 72. The ratio 
of the number of dofs for quarks:bosons:leptons is 72:29:18 with a ~ 5:1 ratio 
between hadrons and leptons. 

STRING EVENT GENERATOR 

SR simulations are implemented as an external process in PYTHIA using the 
UPINIT and UPEVNT subroutines. Simulation parameters include M s , Pr-min 
and cuts on the emitted particles. 
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Chapter 4 



Phenomenology of New Physics 

at the LHC 



In this section we present our analysis of BH, SUSY, graviton and SR events at 
the LHC. Part of the original work in this chapter is based on Refs. [81, 82, 83]. 



4.1 SUSY and BH event analysis 

Our analysis of SUSY and BH events is based on visible energy, missing mo- 
mentum of leptons and hadrons, event shape variables and high-P^ leptons. 
Fig 4.1 shows visible energy, f T and sphericity distribution for the five AT- 
LAS points of Table 1.4, where all SUSY processes except SM Higgs production 
have been implemented. For the purposes of the first part of our analysis, the 
difference between the five points is not significant and any of them can be 
chosen as SUSY benchmark. In the following, we will consider point A. This 
is justified by the fact that point A allows for SUSY Higgs production [149]. 
Since BHs may evaporate into Higgs (see Sect. 1.5 below), a meaningful com- 
parison of SUSY and BH events requires the presence of the Higgs channel 
in both models. Moreover, distinguishability of SUSY and BH events must 
be assessed by minimizing the differences between the two models. Since BH 
events are characterized by up to several TeV of ]f T SUSY points with large 
tf T such as point A, must be considered. 

The parameters for the BH benchmark model are fundamental Planck scale 
M* = 1 TeV, minimum BH mass M m i n = 2 TeV, classical-to-quantum thresh- 
old Qmin = 1 TeV, six extra-dimensions (n = 6) and two-body final decay 
(n p = 2). Particles produced in the initial-radiation phase are removed by 
imposing P? cuts of 5 GeV and 15 GeV for leptons and photons+hadrons, 
respectively [146, 147]. 

In Sect. 4.1.1 we will discuss missing energy and momentum followed by 
an analysis of event shape variables in Sect. 4.1.2. BH events tend to be 
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Figure 4.1: Comparison of visible energy (top left), ]f T (top right) and sphericity 
(bottom) for 10000 events for the five ATLAS points of Table 1.4 (A: filled black 
circles, B: filled red squares, C: filled green triangles, D: open blue circles and E: 
open pink squares). 



more spherical than SUSY events due to the spherical nature of the Hawking 
radiation; this is specially evident for high-mass BHs. The formation of a stable 
BH remnant at the end of the evaporation phase also helps to discriminate 
SUSY and BH events because of the large amount of energy which is carried 
away by the remnant. 

Analysis of isolated dilepton events is another powerful method to distin- 
guish the two models and will be discussed in Sect. 4.1.3. This is due to the 
fact that SUSY dileptons typically originate from a single decay chain whereas 
leptons are rarely emitted by BHs (the hadron-to-lepton ratio is approximately 
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5:1) and are uncorrelated; i.e. they can be emitted at any angle w.r.t. beam 
axis. 

4.1.1 Energy and Momentum 

Figure 4.2 shows visible energy, lf T and Pt of leptons and hadrons & photons 
for 10,000 SUSY and BH benchmark events. The amount of visible energy 
and If T is comparable for the two scenarios, even in the absence of a BH 
remnant. This is due to the presence of invisible channels in both models: 
neutrinos+gravitons for the BH and LSP for SUSY. The flavor of the decay 
products is a better discriminator. SUSY interactions do not produce leptons 
with energy above the TeV since isolated leptons are produced by the decay of 
sparticles with typical energy of less than a few hundred GeV. On the contrary, 
quanta produced in the BH decay are characterized by an average energy 
E ~ M/N, where the multiplicity N is less than 10 for typical BHs at the 
LHC. Since Hawking evaporation does not distinguish leptons from hadrons, 
hard leptons with energy up to several TeV are likely to be produced during 
the BH decay. This suggests that isolated leptons may provide a powerful 
means to discriminate the two models. This is indeed the case, as we shall see 
in Sect. 4.1.3. 

Fig. 4.3 shows how variations in the BH Planck phase affect the observ- 
ables of Fig. 4.2. The plots compare zero (BH remnant)-, two- and four-body 
decays. The BH is expected to have shed its electric and color charges by the 
time the remnant is formed 1 . The BH remnant is thus undetectable and a 
source of f T in addition to neutrinos and gravitons which are emitted during 
the Hawking evaporation phase. This leads to a larger difference in If T between 
SUSY and BH models. The visible Pt in hadrons & photons is sensibly re- 
duced in the presence of a BH remnant. The BH remnant carries away energy 
which otherwise would have been emitted in visible channels (mostly hadrons) 
during the BH decay phase. The leptonic channel is essentially unaffected by 
the presence of a BH remnant since leptons are rarer than hadrons in the BH 
decay phase. Variations in the energy distribution of the leptonic channel are 
thus suppressed compared to the hadronic channel. Changes in the number 
of final Planckian hard quanta do not produce significant differences in the 
distributions; more quanta of lower energy behave statistically like less quanta 
with higher energy. Provided that the BH decays at the end of the Hawk- 
ing phase, it is thus safe to set the number of Planckian quanta to n p = 2 
(or n p — 4), although BHs may decay in different numbers of particles on a 
event-to-event basis. Variations in the classical-to-quantum threshold Q m in 
are also not expected to cause significant differences in the energy/momentum 
distributions. A higher threshold increases the emission in the Planck phase 
while decreasing Hawking radiation. Since these phases differs only in relative 

1 For an alternative scenario see, Ref. [150]. 
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Figure 4.2: Comparison of 10,000 SUSY and BH benchmark events at the LHC. 
Visible energy and If T (top panels) are comparable due to the presence of invisible 
channels in both models. Leptons with large Pt provide instead an effective discrim- 
inator (bottom right panel). 



greybody factors, the effect is too small to be detected. 

Fig. 4.4 shows visible energy, and visible transverse momenta of leptons 
and hadrons & photons for different values of the fundamental Planck scale. 
Higher values of M* lead to more massive BHs, higher multiplicity and more 
energetic quanta. This causes a significant increase in missing and visible 
momenta. If the value of the fundamental Planck scale happens to be large, 
BHs are likely to be detected and easily distinguished from SUSY through 
detection of highly-energetic isolated leptons and hadronic jets. If T of several 
TeV would also be observed. However, the LHC is expected to produce very 
light BHs (see Fig. 1.10) therefore such an analysis may not provide the most 
effective discriminators [146, 147, 151, 152, 153, 154, 155, 156]. The study of 
leptonic final states alleviates this problem. See Sect. 4.1.3. 



4.1.2 Event shape variables 

Event shape variables, such as the sphericity, 2 nd Fox- Wolfram moment and 
thrust, and jet masses defined in Chapter 2 can be used to complement the 
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Figure 4.3: Distribution of visible energy, lf T and transverse momenta of leptons 
and hadrons & photons. SUSY plots are shown as pink open squares. The four plots 
show the effect of different decay modes in the Planck phase of ten- dimensional BHs: 
remnant formation (n p = 0, black filled triangles), two-body decay (n p = 2, red filled 
circles) and four-body decay (n p = 4, green filled squares). The fundamental Planck 
scale is M* = 1 TeV. 



above analysis. Our analysis shows that BH events are more spherical because 
of the nature of Hawking radiation and the production of larger number of 
jets for SUSY decays. Formation of a BH remnant and high values of the 
fundamental scale lead to a significant higher sphericity than SUSY events 
(top panels of Fig. 4.5). The 2 nd Fox- Wolfram moment (middle panels of 
Fig. 4.5) is stable versus changes in the BH Planck phase and provides a good 
SUSY/BH discriminator. BH models with higher M* can be differentiated 
more easily from SUSY events. 

Similar conclusions can be reached by looking at jet masses and number of 
jets. SUSY events generate more and lighter jets than the BH model due to 
copious production of quarks (Fig. 4.6). The difference is again especially sig- 
nificant for high values of M* and in the presence of BH remnants. Absence of 
sub-Q min hard jets could provide strong evidence for BH remnant production. 
(See the suppression of heavy jets below the classical-to-quantum threshold 
Qmin = 2 TeV in the top leftmost panel of Fig. 4.6.) 
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Figure 4.4: Distribution of visible energy, lf T and transverse momenta of leptons 
and hadrons & photons. SUSY plots are shown as pink open squares. The plots 
show the effect of varying the fundamental Planck scale: M± = 1 TeV (black filled 
triangles), M* = 2 TeV (red filled circles) and M* = 3 TeV (green filled squares). 
The ten- dimensional BHs decay in two hard quanta at the end of the evaporation 
phase. 



4.1.3 High-P T leptons 

SUSY events can be categorized by the number of high leptons. Events with 
no isolated leptons can occur when a SUSY decay chain ends in the production 
of jets. Events with a single isolated high P? lepton are difficult to identify 
because of the large SM background from W decays into leptons. Events with 
two isolated high P? leptons are characterized by a lower signal-to-background 
ratio (due to the low branching ratio into leptons) with respect to zero lepton or 
one lepton events. Nevertheless, dileptons provide a clear signature because 
the presence of two isolated leptons can be accurately measured with early 
LHC data. 

Production of gluinos and squarks dominate SUSY production at the LHC. 
The bulk of the decay modes of these sparticles are characterized by the pro- 
duction of the lighter chargino and two neutralinos. Thus, the decay of the sec- 
ond lightest neutralino fo rm s an important part of this analysis. Charginos 
decays are more difficult to analyse as they involve either a missing neutrino or 
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Figure 4.5: Sphericity (top panels), 2 nd Fox-Wolfram moment (middle panels) 
and thrust (bottom panels) for 10,000 BH and SUSY events. The left panels show 
the effect of different Planckian decay modes: BH remnant (black filled triangles), 
two-body decay (red filled circles) and four-body decay (green filled squares). The 
fundamental scale is M+ = 1 TeV and the number of EDs n=6. The right panels 
show the effect of different fundamental scales: M* = 1 TeV (black filled triangles), 
2 TeV (red filled circles) and 3 TeV (green filled squares). The ten- dimensional BHs 
decay in two quanta at the end of the Hawking phase. 



quark jets. Two body decays of neutralinos have large branching ratios. If two 
body decays are not kinematically allowed, then the three body leptonic mode 
dominates [137]. A list of dominant SUSY interactions leading to dileptons at 
the LHC is described below for the LHC points of Table 1.4. 

The relevant decay chains for ATLAS point A are shown in Fig. 4.7. The 
third decay chain allows the separation of leptons from the hadronic back- 
ground [157]. Requiring the leptons to be isolated and have high P T (in excess 
of 15 GeV) ensures that they were indeed produced from the above decay 
chain and not from the initial 2x2 hard interaction. In the rest frame of the 
second lightest neutralino, the dilepton invariant mass is 
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Figure 4.6: Heavy and light jet masses (top and middle panels) and number of 
jets (bottom panels) for 10,000 BH and MSSM events. Symbols are like in previous 
figures. 



The invariant mass is a function of the masses of the particles involved in 
the decay. Since the momentum of the LSP is not constrained, the invariant 
mass distribution shows an edge at ~ 100 GeV. This is because Mu reaches 
a maximum when the LSP is produced at rest. In that case the dilepton 
invariant mass is the difference between the masses of the LSP and x° which 
is ~ 100 GeV for point A. 

Heavy sparticles are produced at points B and C, where m =400 GeV 
and mi/ 2 =400 GeV. The dominant decay mode at these points is x° — > Xih 
because x° decay into sleptons is forbidden due to the high mass of the slepton 
[8, 9, 10]. Dileptons are produced from the decay of heavy charginos. Thus, 
the dilepton invariant mass plot does not have an edge as observed at point 
A. The large value of tan /3 at point C increases production of Z° bosons and 
thus the number of dileptons. 

Point D (m =200 GeV, m 1/2 =100 GeV) is characterized by light SUSY 
particles, thus sparticle production rates are large [149]. The dominant SUSY 
processes, where the percentages in parenthesis are the corresponding BRs, is 
gg with g -> lb (89%) -> $66 (86%) -> x\bbl + l~ (33%) [8, 9, 10]. The above 
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Figure 4.7: The top three SUSY decay chains at point A and their branching ratios. 

interactions imply that sparticle decays at this point are characterized by b 
jets and dileptons. 

At point E (mo=800 GeV and mi/2=200 GeV) gluino production is domi- 
nant because M~ g < Mq. Neutralino decay into sleptons is suppressed because 
of the high mass of the slepton. There is significant Z° production a feature 
common to models with large tan f3. Dominant processes and branching ratios 
at this point are xi A -> ^°X?, 2 (39%) and \ 2 -> Z°xi (32%) [8, 9, 10]. 
Dileptons are produced from x° — ¥ Xil + l~ (6%) [8, 9, 10] and from the decay 
of Z° into leptons. 

Point LM1 [6, 11] (mo=60 GeV and mi/2=250 GeV) is characterized by 
a large value of tan/3, which leads to a lower mass of f and b. This results 
in the decay of the second lightest neutralino into f, producing dileptons via 

the process x° "~ > t 1 ± t t — > t + t~Xi an d X% ~^ r+T ~Xii where the r decays 
in the lepton channel [158, 159, 160]. Moreover, Z° production is enhanced 
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due to higher sbottom production. This is due to the decay chain b\ — > 
X3 4 b, X3 4 — ► Z°Xi 2 [158, 159]. Since the mass of the gluino is larger than the 
squark mass, the dominant process is g — > qq as at point A. The branching 
ratios of the dominant decays are X2 ~^ ?i r=F (46%), Xi ~^ &il (36%) and 
X2 ~^ (H-6%) [11]. The invariant mass distribution of dileptons is expected 
to show an edge due to the process X2 ~~ * IrI- 

The SM dilepton background contribution can be removed by applying 
suitable cuts on the transverse momenta of the leptons [96]: Pti > 15 GeV, 
\r)i\ < 2.5, isolation cut J2i^ti < 7 GeV in a cone of R = 0.2, and ^ T > 200 
GeV. Here Pti is the transverse momentum of the leptons, R = \J At] 2 + A0 2 , 
where <fi and 9 are the azimuthal and polar angles of the lepton w.r.t. beam 
axis. The use of isolated dileptons as a SUSY signature has been previously 
discussed in the literature [149, 161, 162]. Although the production of leptons 
is not as high as the production of colored particles, high energy isolated 
leptons provide a cleaner environment by allowing the removal of the SM 
background. The dominant SUSY interaction for opposite-sign, same-flavor 
(OSSF) dileptons at ATLAS point A is xl ~> l± ~ l ~* l±lT Xi with a 27% 
branching ratio (BR) [157]. The maximum dilepton invariant mass for this 



interaction is 




1/2 

-100 GeV. (4.2) 



Dilepton production in BH events differ greatly from dilepton production 
in the MSSM. There is no single process of dilepton production for BHs. Dilep- 
tons are either produced by the BH directly or by the decay of heavier particles 
such as the Z° boson, ti pairs of a combination of the two. Therefore, the BH 
dilepton invariant mass does not show a sharp cut-off at high energy. Since the 
decay of top quarks into leptons is rare [163], and the BR of Z° into leptons is 
small, r(/ + /^)/r tot ~ 0.034 [1], production of OSSF dileptons is less frequent 
for BHs than SUSY events. Our analysis shows that an OSSF dilepton event 
occurs approximately every 100 BH and 20 SUSY events, with a ~ 1:5 ratio 
of BH-to-SUSY dilepton events at fixed luminosity. 

Fig. 4.8 shows the invariant mass distribution for 2000 SUSY OSSF dilep- 
ton events and 2000 BH OSSF dilepton events. As was expected, the SUSY 
invariant mass distribution shows a sharp edge at ~ 100 GeV. The BH in- 
variant mass distribution is characterized by two peaks, at ~ 90 GeV and a 
smaller peak at 1 TeV, and a tail at high P^- The first peak is due to dilepton 
events from single Z° bosons which are directly emitted by the BH. This is the 
dominant channel of OSSF dilepton production in BH events. The peak at 1 
TeV is due to dileptons emitted at the end of the Hawking phase [146, 147]. 
The BH mass and the number of final hard quanta at the end of the Hawking 
phase have been chosen to be Q m in=^ TeV and n p =2, respectively. Since the 
BH at the end of the Hawking phase is expected to be electrically neutral, 
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Figure 4.8: Invariant mass distribution (in GeV) for 2000 SUSY and BH OSSF 
dilepton events for point A. The SUSY distribution (shaded pink histogram) shows 
the typical endpoint due to the presence of the LSP. The high-Px tail of the BH 
distribution is originated by uncorrelated lepton pairs emitted during the Hawking 
evaporation phase. The final BH decay is in two-quanta. The SM background is 
negligible. 



isolated OSSF dilepton events can occur. This peak is expected to be smeared 
out in a more realistic description of the final BH phase [81]. The high-M^ 
tail of the distribution is originated by pairs of uncorrelated leptons from the 
BH. 

Fig. 4.9 shows the invariant mass distribution for 900 same-sign (SS) SUSY 
and BH events. In SUSY SS dileptons are produced by xf — > H /± // ± ,xf — > 
W ± x < i, and from top quark decay into leptons. SS dileptons from BHs are 
originated either from the BH itself or from the decay of heavier quarks and 
bosons (t,t, W, Z°). In both cases SS dilepton events are expected to be rare 
because of the low BR of W into leptons (~ 11% [1]). SS dileptons are selected 
using the same cuts of the OSSF channel. The high-M^ tail of the BH invariant 
mass distribution can be clearly distinguished from the SS case. Isolation 
requirements reduce the SM background to virtually zero [164, 165]. 

Dilepton events with same sign and/or opposite flavor leptons can also 
be used as discriminators (Fig. 4.10). The "democratic" nature of BH de- 
cay makes events with same/opposite flavor leptons roughly equally probable, 
whereas BRs of SUSY events favor same-flavor dileptons. OS dileptons from 
the decay of X® are of the same flavor. On the other hand, leptonic decays of 
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Figure 4.9: Invariant mass distribution for 900 same-sign dilepton events at point 
A. The SM background is negligible. 

charginos or third generation quarks and squarks would produce equal amounts 
of opposite-sign, opposite flavor (OSOF) dilepton pairs as compared to same 
flavor lepton pairs [14, 142, 143]. 

Table 4.1 shows the BRs of same-/opposite-sign, same-/opposite-flavor iso- 
lated dilepton events for SUSY and BH processes. The dominant channel is 
the OSSF channel for both models. However, SUSY and BH events can be 
easily discriminated by comparing the rate of same-flavor events to the rate of 
opposite-flavor events. 

Table 4.1: BRs of high-P? isolated dileptons for SUSY and BH models. 21,000 and 
100,000 events were simulated in the two cases, respectively, yielding approximately 
1000 dilepton events. OS(SS) stands for opposite (same) sign and OF(SF) denotes 
opposite (same) flavor. 
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The number of isolated, high-P^ leptons can also be used to complement 
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Figure 4.10: Invariant mass distribution for 1100 OSOF dilepton events at point 
A. The SM background is negligible. 

the dilepton analysis (Fig. 4.11). SUSY events are capable of producing up to 
five isolated leptons from the cascade decay of heavy sparticles. Events with 
X2X2 or X1X2 ma y produce four or three isolated leptons, respectively [166]. 
On the contrary, events with three or more isolated leptons are very suppressed 
in BH decays at the LHC energy. Although multilepton events are rare, there 
is very little background and they could be effectively used to distinguish the 
MSSM and the BH model. 

Other effective discriminators can be constructed by looking at dilepton 
events with same sign and/or opposite-flavor leptons. The "democratic" na- 
ture of the BH decay makes all dilepton events roughly equally probable, 
whereas the MSSM favors same-flavor dileptons. Presence of hard opposite- 
flavor leptons is a clear indication of BH decay (Fig. 4.12). Our analysis shows 
that 73% of SUSY dilepton events are OSSF, compared to only 50% in the 
BH model. Conversely, opposite-flavor events are twice more frequent in the 
BH model (40%) compared to the MSSM (21%). 

Dilepton analysis at other LHC points 

Figure 4.13 shows the dilepton invariant mass distributions for points B and 
C. The background cuts remove the SM contribution. Unless otherwise stated 
all dilepton histograms in this section show 2000 OSSF events, 900 SS events 
and 1100 OSOF events for both SUSY and SM. The SUSY invariant mass 
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Figure 4.11: Histogram of the number of events with high-Px leptons for 10,000 
MSSM (pink filled squares) and BH interactions (black open triangles). The number 
of BH events with three isolated leptons is smaller than the number of SUSY events 
by a factor of ~ 20. The probability of producing BH events with four or more 
leptons is virtually zero. 



distribution is represented as the shaded (pink) histogram. 



Table 4.2: BRs of high- Pt isolated dileptons for SUSY points B and C for 50,000 
events. The high value of tan (5 = 1 for point C manifests itself as a greater produc- 
tion of OSSF dileptons. 



High Pt isolated dileptons 
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The dilepton invariant mass distributions for point D (left panel) and E 
(right panel) are shown in Fig. 4.14. From the edge in the OSSF invariant 
mass distribution (top panel of Fig. 4.14) the difference in mass of the two 
lightest neutralinos can be obtained: M^o — M^o ~ 50 GeV. If T calculations 
are generally not performed at point D [14, 142, 143]. This is because lepton 
isolation requirements and the presence of b jets are effective in reducing the 
SM background [8, 9, 10]. The bottom panel of Fig. 4.14 shows the OSOF 
dilepton distribution for ~ 2600 SUSY and BH events. More than half of the 



52 




Figure 4.12: Px scatter plot for ~ iOOO isolated opposite-flavor dilepton events for 
SUSY (pink filled squares) and BHs (black open triangles). BH leptons are harder 
than SUSY leptons and show a larger spread in Pt- 

background contribution is from it events [160]. Simulations show that this 
background is negligible due to the low BR of W bosons into leptons. The 
OSOF dilepton invariant mass distribution at point E is characterized by the 
Z° peak and an undistinct invariant mass edge. The dilepton distribution 
for point LM1 is shown in Fig. 4.15. The edge in the OSSF dilepton invariant 
mass distribution is not as evident as in point A because of the large amount 
of missing momentum which is carried away by the neutrinos in the r decay 
[160]. 

Therefore, in conclusion, it is observed that irrespective of the LHC point 
the dilepton invariant mass always provides effective discrimination between 
SUSY and BHs. 
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Figure 4.13: Invariant mass distributions at point B (left panels) and point C 
(right panels). OSSF (top panels), SS (middle panels) and 0S0F (bottom panels) 
distributions clearly discriminate SUSY and BH events. 
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Figure 4.14: Invariant mass distributions at point D (left panels) and point E (right 
panels). Top Panel: Distribution of 2000 SUSY and BH OSSF dilepton events. At 
point E there is enhanced Z° production because of a high value of tan (3. Middle 
Panel: Event distribution for ~ 1700 SUSY and BH SS dilepton events. The SM 
background is negligible. Bottom Panel: SUSY and BH distribution for ~ 2600 
OSOF events. The major SM background due to it is negligible. 
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Figure 4.15: Top Left Panel: Invariant mass distribution for SUSY and BH OSSF 
dilepton events for the CMS benchmark point LM1. Likewise to point A, the SUSY 
distribution shows the endpoint due to the presence of the LSP. Top Right Panel: 
invariant mass distribution for same sign dileptons. Bottom Panel: OSOF invariant 
mass distribution. A high value of tan f3 enhances the production of OSOF dileptons 
as explained in the text. 
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4.2 ED and BH Event Analysis 



In this section we present our original analysis of BH and graviton events at the 
LHC. Simulation results using energy, momentum and event shape variables 
as discriminators are first presented followed by an analysis of dileptons. 



4.2.1 Energy, Momentum and Event shape variables 
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Figure 4.16: Comparison of 10,000 RS graviton and BH benchmark events. Visible 
energy and ]f T ( top panels ) are comparable due to the presence of invisible channels 
in both models. Bottom panels show visible Pt due to leptons (left panel) and photons 
and hadrons (right panel). 



Figure 4.16 shows visible energy, lf T and Pt of leptons and hadrons & 
photons for 10,000 RS graviton emission and BH benchmark events. The 
amount of visible energy is comparable for the two scenarios, even in the 
absence of a BH remnant. The ^ T is due to the presence of neutrinos for the 
RS graviton model. Graviton decays are also capable of producing isolated 
leptons with high Pp. This explains the peak at 1 TeV in the visible Pt 
distribution of leptons. As with SUSY and BH, isolated leptons may provide 
a discriminating signature between BH and graviton decays (See Sect. 4.2.2.). 
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Figure 4.17: Distribution of visible energy, ]f T and transverse momenta of leptons 
and hadrons & photons. RS graviton plots are shown as pink open squares. The four 
plots show the effect of different decay modes in the Planck phase of ten- dimensional 
BHs: remnant formation (n p = 0, black filled triangles), two-body decay (n p = 2, red 
filled circles) and four-body decay (n p = A, green filled squares). The fundamental 
Planck scale is = 1 TeV. 



Fig. 4.17 compares BH events with zero-, two- and four-body final decay 
to graviton events. The variation of visible energy, and visible transverse 
momenta of leptons and hadrons & photons for different values of the funda- 
mental Planck scale is shown in Fig. 4.18. 

The above analysis is complemented by studying event shape variables. 
Graviton decays are dominated by jets. Formation of a BH remnant and high 
values of the fundamental scale lead to significant higher sphericity than gravi- 
ton decays (top panels of Fig. 4.19). The 2 nd Fox- Wolfram moment (middle 
panels of Fig. 4.19) is stable versus changes in the BH Planck phase and pro- 
vides a good graviton/BH discriminator. BH models with higher M* can be 
differentiated more easily from events with graviton decays. Thrust (bottom 
panels of Fig. 4.19) conveys a similar information. 

Similar conclusions can be reached by looking at jet the masses and at the 
number of jets. Graviton events generate more lighter jets than the BH model 
due to the mass of the graviton being lower than the mass of the BHs (Fig. 
4.20). The difference is specially significant for high values of M* and in the 
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Figure 4.18: Distribution of visible energy, ]f T and transverse momenta of leptons 
and hadrons & photons. RS graviton results are shown as pink open squares. The 
four plots show the effect of varying the fundamental Planck scale: M+ = 1 TeV 
(black filled triangles), M* = 2 TeV (red filled circles) and M* = 3 TeV (green 
filled squares). The ten- dimensional BHs decay in two hard quanta at the end of the 
evaporation phase. 



presence of BH remnants. 



4.2.2 High-P T leptons 

Graviton decay into dileptons provides a cleaner signature. The ATLAS and 
CMS detectors are expected to detect gravitons with mass as high as 3.5 and 
2.08 TeV for c=0.1 and 0.01, respectively [167, 168, 169]. The discovery reach 
is better with electrons than with muons because of a better energy resolution 
of electrons. 

The SM background is removed by imposing the following cuts on leptons to 
remove the SM background: P Tl > 15 GeV, \rji\ < 2.5, isolation cut £\ P^ < 7 
GeV in a cone of R = 0.2 and f T > 100 GeV where the symbols have their 
usual meanings. 

In the RS scenario, the majority of high Pj< leptons originate from the 
decay of the graviton and the Z° boson. The neutral graviton decays into a 
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Figure 4.19: Sphericity (top panels) and 2 nd Fox-Wolfram moment (middle panels) 
and thrust (bottom panels) for 10,000 BH and graviton events (pink open squares). 
The left panels show the effect of different Planckian decay modes: BH remnant 
(black filled triangles), two-body decay (red filled circles) and four-body decay (green 
filled squares). The fundamental scale is M+ = 1 TeV and the number of EDs is six. 
The right panels show the effect of different fundamental scales: M* = 1 TeV (black 
filled triangles), 2 TeV (red filled circles) and 3 TeV (green filled squares). Here, 
the ten- dimensional BHs decay in two quanta at the end of the Hawking phase. 

pair of OS dileptons with a OF-to-SF ratio of ~ 1:14. OSOF lepton pairs are 
originated either from the graviton decay into W bosons or from t quark pairs, 
where the top quark decays into leptons. Production of OSSF dileptons is less 
frequent in the BH model than in the RS graviton model. This is because 
of graviton decays in the Z° channel with the Z° decaying into leptons. Our 
analysis shows a 1:7 ratio of BH-to-ED dilepton events at fixed luminosity. 

Figure 4.21 shows the dilepton invariant mass distribution for RS graviton 
production (shaded pink plot) and the BH model with final two-body decay. 
The RS graviton scenario is characterized by two peaks in the invariant mass 
distribution, at ~ 90 GeV and a second peak at 1 TeV. The first peak is due 
to dilepton events produced from the decay of Z° bosons. The second peak at 
~ 1 TeV is due to graviton decay into a pair of leptons. The tail of the BH 
distribution is due to uncorrelated leptons from the decay phase. Since there 
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Figure 4.20: Heavy and light jet masses (top and middle panels) and number of 
jets (bottom panels) for 10,000 BH and graviton events. Symbols are like in previous 
figures. 



is no similar process in the RS scenario to generate leptons, the invariant mass 
distribution of the graviton does not have a tail at high energies. 

Table 4.3 shows the BR's of 1000 high-Py isolated dileptons for the RS 
graviton model and the BH model. OSSF dileptons are dominant in both 
cases. However, due to the flavor conserving nature of the graviton decay, the 
number of SS dileptons produced in the RS scenario is virtually zero. This is 
not the case for BH events, where SS dileptons account for 18% of the total 
dileptons produced. Therefore, the SS dilepton invariant mass would be a 
promising discriminating signature for RS graviton events and BH events. 

It is worthwhile to study the decay of the Z' boson which is predicted by 
many extensions of the SM [1, 170, 171]. The Z' boson is a heavy, (mass > 500 
GeV) [1], colorless, neutral, spin 1 particle. As a result, Z' decay into leptons 
could compete with graviton decay in the same channel 2 . This is especially 
significant for the LHC which could probe Z' masses in the TeV range. The 
observation of a TeV resonance in the dilepton and diphoton channels could 

2 Howcvcr, the decay of Z' into photons is suppressed. 
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Figure 4.21: Invariant mass distribution of 1100 OSOF dilepton events. Most 
high-Px leptons in the RS model come from the decay of the graviton and Z° bosons. 
This is indicated by the two peaks in the invariant mass distribution. As discussed 
before in Sect. J^.l, the BH distribution is also characterized by two peaks. The 
discriminating signature is the high-P? tail in the BH invariant mass distribution. 



provide a strong evidence of a RS graviton, distinguishing it from a Z' boson. 
An additional discriminating property between the graviton and the Z' is the 
angular distribution of the emitted particles, which is dependent on the spin 
of the resonance [171]. The angular distribution is strongly dependent on the 
production mechanism. The angular distributions of the G* and Z' are shown 
in Table 4.4 [167] and the corresponding distributions are shown in Fig. 4.22. 
Gluon fusion predominates at the LHC. For graviton production, the quark- 
antiquark fusion contribution tends to flatten the angular distribution[168]. 
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Table 4.3: BRs of high-Pr isolated dileptons for the RS and the BH models. Re- 
spectively 2.8 x 10 4 and 10 6 events were simulated in the two cases, yielding approx- 
imately 1000 dilepton events. 



High Pt isolated dileptons 


Graviton 


% 


BH 


% 


OSSF 


987 


99 


523 


50 


OSOF 


13 


1 


341 


32 


ss 








190 


18 



Table 4.4: Angular distributions of Z' boson and graviton decay into dielectrons. 
8 is the angle between the beam direction and the electron in the rest frame of the 
resonance. 



Process 


Distribution 


gg — > G* — > e + e~ 


i - cos 4 e 


qq — > G* — > e + e~ 


1 - 3 cos 2 + 4 cos 4 


99, 99 -> Z' ^> e + e- 


1 + cos 2 
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4.3 SR and SM Event Analysis 



Events with high-P^ photons are a powerful discriminator of SR events from 
the SM background. SR events are characterized by high values of visible 
energy and P T as the photon and the jet originate directly from the 2 x 2 in- 
teraction. Isolated photons provide a further means to extract string signals: 
being directly produced from the SR, they are harder than SM photons. Anal- 
ysis of event shape variables enhance the effectiveness of the above method. 
String simulations use the following benchmark: M s —1 TeV and Pr m j n =50 
GeV. Our analysis uses ~ 10000 and 200 string events for M s =l and 2 TeV 
produced from a sample run of 10 7 events at an integrated LHC luminosity of 
1 fb~ l , respectively. 

Figure 4.23 shows the visible energy (top panel) and the Pt of hadrons & 
photons (bottom panel) for 10 million string+SM and SM only events. The 
visible energy and the Pt are produced by the hard photons and the jets of 
the string decay. Their distributions are characterized by a long tail. The 
observation of events with Pt (visible energy) in excess of 3 (6) TeV would 
provide strong evidence of the existence of a SR. 

In our analysis of the photon momentum, we select the hadronic jets ac- 
cording to the following criteria. The detector is assumed to have an absolute 
value of pseudorapidity of 2.6. This ensures that the jets are originated in the 
hard 2x2 scattering rather than in multiple interactions or from the beam 
remnants. The contribution of jets which do not originate in the hard scat- 
tering are minimized by fixing the minimum transverse energy of all particles 
comprising the jet (EjPrJ to 40 GeV [172]. The particles of the jet must be 
within a cone of R = (Ar/ 2 + A0 2 ) = 0.2 from the jet initiator, where 6 
and 4> are the azimuthal and polar angles of the particle w.r.t. the beam axis, 
respectively. Following Ref. [172], the cuts on the photon jets are Pr 7 > 80 
GeV, r] < 2.6 and an isolation cut £ n Pr < 7 GeV in a cone of R = 0.4. 

The top panel of Fig. 4.24 shows the distribution of the highest Pr of 
isolated photons for string+SM and SM only events. Single photons from SR 
are expected to have a higher Pt than SM photons because they are the direct 
products of the string decay. The main source of background for direct photons 
are jet fluctuations and photons originating from the initial and final state 
radiation [172]. In the former jet consists of a few particles including 

high-Pr 7To's which decay into a pair of photons with a ~ 99% BR [172]. Due 
to the high boost, the photons have a relatively small angular separation and 
therefore "fake" a single photon in the electromagnetic calorimeter. The rate 
of this process is 1 out of ~ 10 3 to 10 4 events [172]. Other sources of fake 
photons are H — > 77 [173] or processes from other exotic phenomena, e.g. 
SUSY [174] or LEDs [88, 99, 100]. Isolation cuts on the photon can effectively 
reduce the number of fake photons. 

The invariant mass plot of the two highest Pt jets for each event is shown 
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Figure 4.23: Top panel: Visible energy distribution for string+SM and SM only 
events. The result for string events is shown by the solid red histogram (M s =l 
TeV) and the dashed green histogram (M s =2 TeV). String events can be identified 
from the high-Ex tail for M s =l TeV. Bottom panel: Distribution of visible Pt for 
^y+hadrons. The high-Px tail is a strong indicator of the presence of SR. 



in the bottom panel of Fig. 4.24. Due to the nature of the interaction under 
consideration, the bulk of the events for both string+SM and SM only events 
are comprised of dijets. These were selected using the same cuts as before. As 
is expected, the SM invariant mass distribution is negligible beyond a certain 
threshold (~ 4 TeV). This is due to the production of direct soft photons and 
jets from the SM interaction. The string+SM distribution is characterized by a 
long tail of energy up to several TeV (three times more than the SM). This tail 
is originated from the decay of string resonances into hard jets and photons. 
Therefore, the measure of a large invariant mass would provide strong evidence 
of a string mediated interaction. 

Histograms for different event shape variables are shown in Fig. 4.25. String 
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Figure 4.24: Top panel: Distribution of the highest JFV for each event. 7 's with 
high Pt created in the string decay are the source of the long tail. Bottom panel: 
Dijet invariant mass distribution. String decays may result is large invariant mass. 



events tend to be more spherical than the background: string+SM decays 
generally produce a distribution of high Pt jets at slightly higher values than 
the SM background. The jets originate from the decay of SRs into photons 
and hadrons. The SM generates less heavier jets than its string counterpart 
as shown in the middle and bottom panels of Fig. 4.25. 

Before ending this section, let us briefly discuss the discovery reach of 
SRs at the LHC Fig. 4.26 shows the signal-to-noise (SNR) ratio for two val- 
ues of the integrated luminosity as a function of M* when Pr min =50 GeV. 
The instantaneous luminosity is assumed to be the LHC design luminosity of 
lO" 34 ^" 2 ^^ 1 corresponding to an integrated luminosity of ~ 1300 (300) fb" 1 
in 4 (1) years. The SNR ratio for an integrated luminosity of 300 fb~ l is ~ 
10 times less than the value obtained at 1300 fb^ 1 . In an optimistic scenario, 
SRs at the TeV scale could be discovered after a successful year of running the 
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Figure 4.25: Histograms of event shape variables for 10 million string+SM and 
SM only events. String events are shown in solid red (M s =l TeV) and dashed green 
(M s =2 TeV). SM events are shown in dotted blue. String events have on the average 
a slightly higher sphericity than SM events due to the slight increase in the number 
of jets (top panel). Similar conclusions are reached from the heavy and low jet mass 
distributions (middle and bottom panels, respectively). 



LHC, assuming a signal-to noise ratio of 10. 
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Figure 4.26: SNR for a integrated luminosity of 1300 fb~ l (solid red line) and 300 
fb^ 1 (dashed green line), corresponding to 4 year and 1 year at an instantaneous 
luminosity of 10 34 cm _2 s _1 , respectively. 

4.4 SR and BH Event Analysis 

A powerful way of discriminating BH and string events would be searching for 
a Z° mass peak in the invariant mass of high-P^ leptons. Z° production is 
highly suppressed in case of string events [129]. On the contrary BH decay 
is characterized by the production of a variety of particles with high P?. A 
rough counting of the number of degrees of freedom of these particles shows 
that the estimated rate of hadron-to-lepton production is 5:1 and the rate of 
Z° and 7 production is comparable (~ 2% to 3%) with the Z° bosons decaying 
into opposite-sign leptons with a 3.4% BR. The invariant mass distribution of 
BH events peaks at ~ 92 GeV, confirming the production of a Z° boson [81]. 
Therefore, the presence of a peak at ~ 92 GeV in the invariant mass of leptons 
would effectively rule out formation of SRs in favor of BH production. Due to 
the nature of the Hawking radiation, BH events are also expected to have on 
average higher sphericity and a larger number of jets than string events. This 
is because in the pp — > 7 + jet channel string events are dominated by dijets, 
whereas BH decay produces multijets. 
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Conclusions 



In this thesis we studied different scenarios for physics beyond the SM: SUSY, 
EDs and SRs. 

Signatures of SUSY and TeV-scale BH events at the LHC have been dis- 
cussed in Sect. 4.1. A combined analysis of event-shape variables and dilepton 
events has shown that it is possible to distinguish BHs and SUSY at the AT- 
LAS and CMS detectors. Event shape variables alone cannot unequivocally 
discriminate between SUSY and BHs. For example, the thrust distribution for 
BH events with no remnant production is indistinguishable from the thrust dis- 
tribution for SUSY events. However event shape variables may prove useful 
when combined with the analysis of the leptonic channel. Isolated dileptons 
could provide the "smoking gun" for detecting BHs at the LHC. The BH dilep- 
ton invariant mass shows a tail at high energy which is absent in the SM or 
MSSM. This analysis can be further strengthened by looking at the number 
and flavor of isolated leptons. 

The possibility of discriminating BH and graviton events was discussed in 
Sect. 4.2. Our simulations showed that by analyzing the dilepton channel, BH 
and graviton events can be clearly discriminated at the LHC. As with SUSY 
and BHs, event shape variables provide additional means of discrimination. 

Creation of SRs at the LHC was investigated in Sect. 4.3. If the string 
scale is ~1 TeV, our analysis has shown that SRs can be detected. With 
the LHC operating at a luminosity of 2 x 10 33 cm" 2 s~ 1 SRs could be detected 
within a year of the start of the LHC. String events show higher sphericity 
and higher visible energy than the SM background. These quantities provide 
an effective means of detection when combined with the measure of the Pt 
of isolated photons and the dijet invariant mass. Since the final products of 
the SR are directly produced from the string decay, the dijet invariant mass is 
characterized by a tail at high energies which is absent in the SM. A powerful 
way of discriminating between BH and string events would be searching for a 
Z° mass peak in the invariant mass of high-P^ leptons. SR decay into leptons 
and photons may be another effective method for distinguishing SRs and BHs. 

The investigations in this thesis provide a starting point for a 
quantitative study of discrimination of BHs from SUSY, graviton events and 
SRs at the LHC. However, many open questions remain unanswered. For 
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example, our work has not dealt in depth with sources of background at the 
LHC. The actual sensitivity of the method for specific experiments (CMS or 
ATLAS) is also not discussed here. Moreover, the kinematical cuts 
implemented in Chapter 4 although commonly used in the literature, lack a 
study at the generator level of their efficiency. Finally, it would also be 
worthwhile to study the detector response with CMSSW, the software 
framework which is used for simulation and analysis of high energy collisions 
at the CMS detector. 
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Appendix A: Units 



In this thesis we use 3 different system of units: high energy, SI and natural. 
Table 1 shows the value of the fundamental constants in the 3 systems of units. 
Conversion factors between these 3 systems are summarised in Table 2. 



Table 1: Fundamental constants. 



Quantity 


High energy units 


SI units 


Natural units 


h 


6.588XKT 28 GeV s 


1.055xl0~ 34 J s 


1 


c 


2.998 xlO 23 fm/s 


2.998x10 s m/s 


1 


G4 


6.707xl(r 39 GeV 


6.674xKr n m 3 kg~V 2 





Table 2: Conversion factors. 



Quantity 


High energy units 


SI units 


Natural units 


Length 


1.97xl0 -20 fm 


1.97xl0" 35 m 


1 M PL 


Mass 


1.22xl0 19 GeV 


1.78xl(T 8 kg 


1 M PL 


Energy 


1.22xl0 19 GeV 


1.95xl0 9 J 


1 M PL 


Time 


8.2xl(T 2U GeV" 1 


6.58 xKT 44 s 


1 Mp I 


Temperature 


1.22xl0 19 GeV 


1.42xl0 32 K 


1 M PL 
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Appendix B: Acronyms 



ADD 


A 1 *TT J T\* 1 1 ' 

ArkamHamed-Dimopoulos-Dvah 


BD 


Black Disk 


T T 

BH 


1 1 XT 1 

BlackHole 


BK 


Branching Ratio 


CERN 


European Organization tor JNuclear Research 


CL 


Confidence Level 


CM 


Center-of-Mass 


CP 


Charge-Parity 


dor 


degrees ot treedom 


ED 


Extra Dimension 


GUT 


iTT ■ f* i mi 

Grand Unmed Iheory 


KK 


Kaluza-Klein 


LED 


Large Extra Dimension 


LEP 


Large Electron-Positron 


LHC 


Large Hadron Collider 


t on 

LbP 


Lightest Supersymmetnc rarticle 


MC 


A/Tontp C^rlo 


MSSM 


Minimal Supersymmetric extension of the SM 


mSUGRA 


minimum SUperGRAvity 


OSOF 


Opposite Sign Opposite Flavor 


OSSF 


Opposite Sign Same Flavor 


PDF 


Parton Distribution Function 


QCD 


Quantum ChromoDynamics 


QED 


Quantum ElectroDynamics 


RS 


Randall- Sundrum 


SM 


Standard Model 


SR 


String Resonance 


SSOF 


Same Sign Opposite Flavor 


SSSF 


Same Sign Same Flavor 


ST 


String Theory 


SUSY 


SUperSYmmetry 


TS 


Trapped Surface 



N6 



Appendix C: Symbols 



rribare Bare Higgs mass 

rriHiggs Physical Higgs mass 

A (Sect. 1.2) Higgs mass cutoff 

X and \x (Sect. 1.2) Momentum transfer scale 

hi Numerical prefactors 

mo common mass of scalar particles at Mqut 

mi common mass of gauginos and Higgsino at M GUT 

2 

A trilinear coupling at Mqut 

tan f3 ratio of vacum expectation values of two Higgs fields 

fi sign of the Higgsino mass parameter 

L(B) Lepton (baryon) number 

s Particle spin 

n Number of extra-spatial dimensions 

R Size of ED in ADD model 

ip Scalar field 

m^p Mass of the scalar field 

Xi 3 spatial coordinates 

G4 4-dimensional Newton's constant 

G n+ 4 (n+4)-dimensional Newton's constant 

Sd Surface area of unit sphere in d dimensions 

F 4 Force between two test masses in 4 dimensions 

F n+4 Force between two test masses in (n+4) dimensions 

D Total number of dimensions 

4> g {i) 4-dimensional gravitational potential 

^g(n+4) (n+4)-dimensional gravitational potential 

p m Mass density 

V n Volume of ED in ADD model 

x n Bessel function roots 

A (Sect. 1.4.3) Strength of graviton-matter coupling in RS model 

a cross section 

Mpl Planck mass 

M BH mass 

s Parton square CM energy 

s Particle square CM energy 
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Q (Sect. 1.5) 


Momentum transfer scale 


T H 


Hawking temperature 


T 


BH lifetime 


i 

a 


Slope parameter 


Is 


String length scale 


M s 


String mass scale 


9s 


String coupling 


9 W 


Weinberg angle 


Pz 


Momentum along z axis 


L 


Luminosity 


f 


Proton bunch crossing frequency 


n b 


Number of protons in each bunch 


A 


Bunch area 


S 


Sphericity 


T 


Thrust 


R2 


Second Fox- Wolfram moment 




Visible energy 


Pt 


Transverse momenta 


p 

1 Train 


Minimum transverse momenta 


Pti 


Transverse momenta of leptons 


Pt 


Transverse momenta of photons 


¥t 


Missing transvese momenta 


$T 


Missing transvese energy 


Tj 


Pseudorapidity 
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